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TASK  TITLE;  Analysis  of  Impurity  Band  Formation  and  Related  Effects 
Near  the  P/N  Junction  in  Compensated  GaAs  and  Si. 


TASK  OBJECTIVE;  Utilizing  the  bandgap  narrowing  (BGN)  model  developed 
earlier  by  the  author,  the  objective  is  to  apply  the  model  in 
studying  concentration  dependent  changes  near  the  P/N  junction 
of  semiconductor  parameters  that  have  a  direct  effect  on 
current-voltage  characteristics. 

TASK  DESCRIPTION;  Based  on  statistical  variations  in  neighboring 
impurity  separations,  free-carrier  screening  of  impurity 
interaction  energies,  and  statistical  fluctuations  in  the 
local  impurity  concentration,  the  BGN  model  is  to  be  used  in 
studying  variation  across  the  P/N  junction  of  several  important 
characteristics  including  bandgap  narrowing  and  edge  tailing, 
impurity  band  tailing  and  merging  with  free  carrier  bands, 
Fermi-level,  n-p  type  reversal  due  to  fluctuations  in  local 
impurity  density,  which  may  influence  recombination  rate. 
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program  was  sponsored  by  the  Air  Force  Wright  Aeronautical  Laboratories 
Wright-Patterson  Air  Force  Base,  Ohio.  The  Air  Force  Project  Engineers 
were  W.  P.  Rahilly  and  Joseph  Wise  (AFWAL/POOC) . 
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1 .  Introduction 


Fabrication  techniques,  used  in  the  production  of  photovoltaic  devices 
such  as  solar  cells,  usually  lead  to  high  donor  and  acceptor  concentrations, 
especially  in  compensated  materials,  sufficient  to  introduce  substantial 
interactions  between  neighboring  impurities.  Because  of  these  interactions, 
and  also  because  of  the  electric  screening  effect  of  free  charge  carriers, 
the  donor  and  acceptor  energy  levels  are  shifted  with  respect  to  the  edge 
of  the  conduction  or  valence  band,  forming  impurity  bands  and  causing  changes 
in  the  band  gap.  In  addition  to  the  obvious  variations  with  temperature  and 
impurity  concentration,  these  effects  also  vary  from  region  to  region  within 
the  semiconductor  material,  due  to  the  statistically  random  nature  of  dopant 
implantation,  and  vary  with  usage  due  to  attendant  changes  in  charge  carrier 
density  near  the  P/N  junction. 

A  comprehensive  analytical  treatment  of  these  effects  is  not  feasible, 
primarily  because  the  energy  interactions  between  impurity  ions  and  localized 
charge  carriers  constitute  a  many-body  problem  in  quantum  theory,  but  also 
because  of  the  random  statistical  nature  of  the  positions  of  impurity  ions  in 
the  host  lattice.  Many  investigators  have  been  active  [1  -  24]  in  various 
aspects  of  the  overall  problem,  however,  no  treatment  has  been  found  which 
energy  levels  of  interacting  impurities  are  calculated  employing  a  screened 
potential  energy.  Some  studies  [16,  33]  have  been  made  with  density  matrix 
theory  and  potential  energy  screening  leading  directly  to  density-of-state 
data  but  these  do  not  provide  detailed  information  about  individual  energy 
levels  needed  for  developing  insight  into  the  behavior  of  the  semiconductor 
materials  in  question.  In  most  studies,  one  or  more  of  the  above  effects 
have  been  treated  only  qualitatively  or  not  at  all. 

A  compromise  approach  is  used  in  the  present  study.  In  order  to  make 
the  analysis  feasible,  only  nearest-neighbor  impurity  ions  are  considered, 
but  the  energy  levels  are  calculated  with  a  screened  potential  interaction. 


The  resulting  levels  are  then  used,  along  with  probability  distributions  to 
account  for  the  random  positions  of  the  impurity  ions,  to  calculate  densities 
of  state  for  each  of  the  levels  considered.  In  this  way,  the  most  important 
contributing  influences,  interactions  between  neighboring  impurities,  charge 
carrier  potential  screening,  and  random  distances  between  impurities,  are  all 
treated  quantitatively,  thus  providing  a  procedure  for  following  changes  in 
device  characteristics  due  to  changes  in  temperature,  impurity  concentration, 
compensation,  and  charge  carrier  densities. 

In  addition  to  the  physical  characteristics  of  the  host  material,  the 
input  parameters  needed  are  temperature,  donor  concentration,  and  acceptor 
concentration.  Before  the  energy  levels  can  be  calculated,  the  Fermi  level, 
charge  carrier  density  and  screening  factor  are  needed  and  must  be  estimated 
to  provide  this  input.  After  the  energy  levels  are  then  determined,  they 
are  used  to  recalculate  the  estimated  values  above  and  the  energy  levels 
recalculated.  This  procedure  is  repeated  until  a  self-consistent  set  of 
values  is  obtained.  The  energy  levels  calculated  are  those  for  a  pair  of 
impurity  ions  separated  by  a  distance,  R,  similar  to  the  treatment  of  a 
diatomic  molecule  except  that  R  is  constant  for  a  given  pair,  and  also  the 
localized  carriers,  bound  to  the  ions,  are  assumed  to  be  under  the  influence 
of  a  screened  electric  potential  instead  of  the  usual  Coulomb  potential. 
Levels  are  calculated  for  each  of  eight  different  systems  involving  pairs, 
donor-donor,  donor-acceptor,  and  acceptor-acceptor  with  either  one  or  two 
localized  charge  carriers.  The  levels  are  obtained  for  a  range  of  values 
for  the  distance,  R,  to  cover  all  possibilities  needed  for  the  statistically 
random  distances  found  between  nearest  neighbors.  The  probability  of 
finding  a  given  R  is  then  calculated  as  a  function  of  R  and  used,  along  with 
the  energy  levels  to  obtain  the  probability  of  finding  a  given  energy  level 
as  a  function  of  energy. 

These  calculations  are  performed  for  a  set  of  temperatures  over  values 
of  interest  and  for  a  range  of  values  of  donor  concentration  as  well  as  for 
acceptor  concentration.  The  resulting  information  is  presented  in  the  form  of 
graphic  displays  found  in  Appendix  A. 
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2.  Physical  Model 

Eight  different  physical  systems  are  involved  in  the  calculation  of 
energy  levels.  Four  of  these  form  donor  levels  labeled  E^ ,  E 2,  E^,  4  E^, 
as  follows:  E1  is  the  groundstate  energy  of  a  single  electron  bound  to 
a  pair  of  donor  ions,  abbreviated  (DD)-  ;  Ez  is  the  additional  energy  of 
the  second  electron  bound  to  a  donor  pair  with  the  first  electron  present, 
abbreviated  <DD->-  ;  E^  is  the  groundstate  energy  of  an  electron  bound 
to  a  donor-acceptor  pair  with  no  hole  present,  abbreviated  (DA)-  ; 

E^  is  the  energy  of  an  electron  bound  to  a  donor-acceptor  pair  with  a 
hole  present,  abbreviated  (DA+)-  .  The  other  four  systems  form  acceptor 
levels  labeled  ,  H 2,  Hg,  &  H^,  abbreviated  (AA)+  ,  (AA+)+  ,  (AD)+,  & 
(AD-)+  ,  respectively.  For  the  pair  separation,  R,  large  compared  to  the 
effective  Bohr  radius,  the  above  systems  reduce  to  isolated  systems,  for 
example,  E^^  becomes  the  energy  of  one  donor  atom  (electron  present)  and 
one  donor  ion,  a  total  of  one  effective  Rydberg  for  an  electron.  The 
binding  energy  of  an  isolated  donor  electron  is  not  necessarily  the  same  as 
that  of  a  hole  to  an  isolated  acceptor.  In  this  study,  the  isolated  donor 
is  taken  as  the  unit  of  energy,  one  Rydberg,  and  the  isolated  acceptor  is 
assumed  to  have  a  different  groundstate  hole  energy  in  order  to  display 
in  the  results  the  consequences  of  an  inequality.  This  introduces  a 
corresponding  inequality  in  the  effective  Bohr  radii  with  that  of  the 
acceptor  being  different  in  inverse  proportion  to  the  energies.  The  ratio 
of  these  energies  is  an  input  parameter  of  the  computer  program. 

The  donor  energy  level  calculations  follow  the  procedures  used  in  the 
treatment  of  the  hydrogen  molecule  ion,  ,  in  the  cases  of  levels,  and 
the  hydrogen  molecule,  H^,  in  determining  E^  for  two  electrons,  from  which, 
E 2  *  “  E, ,  is  found.  A  similar  procedure  is  used  for  the  acceptor 

levels  and  in  both  cases  a  variation  procedure  is  used  to  minimize  the 
energy  obtained  with  s-state  basis  wavefunctions,  the  details  are  contained 

3 


CLLU I'vlvlv 


■-V- 


V  V 


in  Appendices  B  and  C.  In  the  case  of  the  donor-acceptor  pair,  a  modified 
but  similar  calculation  is  developed  to  account  for  one  positively  and  one 
negatively  charged  ion  in  the  pair.  Two  variation  parameters  are  used  in 
order  to  obtain  better  approximations  to  the  energies  which  are  minimized 
with  respect  to  a  variation  of  both  parameters.  The  details  appear  in 
Appendix  B  in  which  and  E^  are  also  called  and  A  similar  method 

is  used  to  obtain  values  for  H-  and  H.. 

3  4 

In  the  hydrogen-like  calculations,  certain  modifications  are  needed 
in  order  to  properly  represent  the  change  from  a  hydrogen  molecule  in  space 
to  a  hydrogen-like  system  in  the  crystal  lattice: 

<1>  The  screened  potential,  ±2  exp(-f3r)/r  ,  is  used  instead  of  the  Coulomb 
potential,  ±2/r  ,  where  (3  is  the  screening  factor  (also  labeled  B).  The 
units  are  effective  Bohr  radii  and  effective  Rydbergs  for  the  isolated 
donor  electron.  The  potential  is  assumed  to  be  isotropic  in  order  to 
keep  the  calculation  manageable. 

(2)  The  dielectric  constant,  the  electron  effective  mass,  and  the  hole 
effective  mass  are  assumed  to  be  isotropic  and  are  used  to  modify  the 
Bohr  radius  and  the  Rydberg  to  their  effective  values  in  the  lattice. 

(3)  The  potential  energy  term,  ±2/R  ,  between  the  impurity  ions  is  not 
included  in  the  Hamiltonian  since  the  ion  separation,  R,  is  not  a 
dynamic  variable  (in  contrast  with  the  hydrogen  molecule).  The  value 
of  R  may  vary  from  one  ion  pair  to  another,  because  of  the  random 
locations  of  impurities,  but,  in  a  given  pair,  the  value  of  R  is  fixed. 
For  each  of  the  values  of  temperature,  donor  and  acceptor  concentration 

considered,  the  eight  energy  levels  are  calculated  for  values  of  R  covering 
a  range  up  to  eight  Bohr  radii  in  32  steps.  The  first  such  calculation  is 
based  on  a  value  of  the  screening  factor,  (3,  obtained  from  a  simple  estimate 
of  the  Fermi  energy  and  the  free  charge  carrier  density  using  the  Thomas- 
Fermi  approximation  for  (3.  The  resulting  energy  levels  are  then  used  to 
obtain  an  improved  value  for  (3  based  on  the  distribution  of  energy  levels 
and  their  probability  of  occurrence  (see  Appendix  B  for  details).  The  energy 
levels  are  then  recalculated  and  this  procedure  is  repeated  until  no  futher 


changes  occur  in  the  value  of  13,  thus  assuring  self-consistent  results. 

In  the  recalculation  of  the  Fermi  energy,  E_,  and  the  screening  factor, 

4-  * 

(3,  the  charge  balance  equation,  n+N  =  p  +  N  ,  is  used  thus  insuring  the 

results  to  be  valid  for  equilibrium  conditions  in  bulk  material.  Also  used 

in  this  step  is  the  probability  distribution  function  which  describes  the 

relative  frequency  of  occurrence  of  different  values  of  the  distance,  R, 

between  neighboring  impurity  ions.  The  relation  follows  Poisson  statistics, 

u 

i.e.,  Pk<cv)  =  <(cv)  / k!)  exp(-cv)  ,  which  in  effect  states  that,  if  the 

average  number  of  impurity  ions  expected  in  the  volume,  v,  is  given  by  cv, 

then  the  probability  of  finding,  instead,  the  number,  k,  is  given  by  the 

expression  above.  The  validity  of  this  formulation  depends  on  only  the 

facts  that  the  position  of  an  ion  in  the  lattice  is  not  dependent  on  the 

positions  of  the  other  ions  and  the  number  of  impurity  ions  is  much  smaller 

than  the  number  of  lattice  sites.  The  distribution  of  nearest  neighbor 

impurity  ion  separations,  R,  can  be  deduced  from  this  expression  to  be  given 

2  3 

by,  P(R,c')dR  =  4xrc'R  dR  exp(-4nc'R  /3>  ,  where  R  is  in  units  of  Bohr  radii 
and  c'  is  the  impurity  concentration  in  the  corresponding  units,  namely, 

-3 

(Bohr  radii) 

In  compensated  semiconductors  with  a  donor  concentration  of  in  the 
bulk  material,  and  an  acceptor  concentration  of  diffused  into  the  device 
from  a  surface  value  of  Cg  ,  the  concentration,  c',  in  the  probability  above 
is  derived  from  the  total  impurity  concentration,  C.  =  C  +  C.  .  The  chance 

v  X  D  2 

of  finding  D-D,  D-A,  &  A-A  pairs  at  distance  R  is  then  proportional  to,  q  , 
2q(l-q),  and  (l-q)(l-q),  respectively,  where  q  =  C^/C^.  The  products  of  the 
probability  above  and  these  factors  then  provide  the  necessary  information 
for  weighting  the  different  energy  levels  in  the  process  of  determining  the 
contributions  to  the  charge  balance  equation,  above,  and  in  generating  the 
desired  density-of-states  data.  More  details  of  this  process  can  be  found 
in  Appendix  B. 

Five  temperatures  were  chosen  for  consideration  in  the  computation, 
namely,  100  K,  200  K,  300  K,  400  K,  and  500  K.  For  each  of  these  values, 
five  bulk  (background)  donor  concentrations  were  selected.  Their  values 


are  Log  C.  =  15.5,  16,  16.5,  17,  and  17.5  ,  or  C.  =  3.16x10*^,  10*®, 

16  C  17  17  -3  “ 

3.16x10  ,  10  ,  and  3.16x10  cm  .  For  each  of  these  cases,  nine  values 

of  the  variable  acceptor  concentration,  C  ,  were  selected  to  cover  a  range 

of  positions  in  the  vicinity  of  the  P/N  junction.  Figures  1  (a)  and  1  (b) 

show  typical  concentration  profiles  across  a  P/N  junction  in  which  the 

Log  is  assumed  to  follow  a  parabolic  decrease  in  magnitude  starting  at 

a  value  Log  at  the  device  surface,  and  the  junction  is  assumed  to  be 

located  at  a  distance  of  0.3  or  0.4  microns  from  the  surface. 

The  values  of  C  and  C.  are  shown  as  dotted  lines,  C.  is  shown  as  a 
x  d  t 

dashed  line  and  the  net  concentration,  C  =  IC^  -  C^l  ,  is  shown  as  a  solid 
line.  The  electrical  characteristics  of  the  device  are  closely  related  to 
the  net  concentration,  C,  however,  those  effects  which  result  from  closely 
spaced  impurity  ions  also  depend  strongly  on  the  total  concentration,  C^.. 
Nine  numbers,  0-8,  serve  to  mark  the  points  on  the  curve,  C,  which  follow 
from  the  selected  values  of  mentioned  above.  Actually,  the  values  of  C 
were  chosen  to  be  equal  to  (0  and  8)  or  such  that  Log  C  was  less  than 
Log  Cb  by  steps  of  0.5  ,  1.0  ,  1.5  ,  and  2.0  ,  on  each  side  of  the  junction 
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The  numerical  results  obtained  by  the  use  of  this  model  are  presented 

in  Appendix  A  in  the  form  of  graphic  displays.  Each  figure  consists  of  five 

displays,  (a)  Impurity  Spacing  Probability  Distributions,  <b>  Donor  Energy 

Levels,  (c)  Acceptor  Energy  Levels,  (d)  Densities  of  States  for  Conduction 

Band  and  Donor  Levels,  and  <e)  Densities  of  States  for  Acceptor  Levels  and 

Valence  Band.  Each  Fig. (page)  is  based  on  input  data  for  a  single  selection 

of  temperature,  T,  donor,  ,  and  acceptor,  C^  ,  concentration.  There  are 

five  temperatures  treated,  100  K,  200  K,  300  K,  400  K,  and  SOO  K.  The  first 

digit  of  the  figure  number  is  the  first  digit  of  the  temperature. 

For  each  temperature,  there  are  five  selections  of  background  (or  bulk) 

dopant  level  treated  such  that  Log  C.  =  IS. 5,  IE,  IE. 5,  17,  and  17.5,  where 
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C.  is  in  units  of  cm  C.  and  C  are  also  labeled  Cl  and  C2.  The  second 
b  b  x 

digit  of  the  figure  number  has  values  0,  1,  2,  3,  and  4,  corresponding  to 

the  values  of  Log  Cfa  above,  respectively. 

For  each  of  these,  nine  values  of  the  infused  dopant,  C^,  are  chosen  to 

form  a  profile  of  values  across  the  P/H  junction  (see  points  labled  0-8 

shown  in  Fig.  1,  Typical  Concentration  Profiles).  These  labels  correspond 

to  the  last  digit  in  the  figure  numbers.  The  points  0,  1,  2,  and  3  are 

chosen  so  that  Log  C  is  less  than  Log  C^  by  values  of  0.0,  0.5,  1.0  and  1.5, 

where  C  =  C  -  C.  .  The  point  labeled  4  is  at  C  *  C.  -  C  =  1010. 

x  b  b  x 

Points  5,  £  and  7  provide  values  of  Log  C  less  than  Log  Cfc  by  1.5,  1.0 

and  0.5,  respectively,  with  C  *  C.  -  C  ,  and  point  8  sets  Log  C  almost 

0  -6  * 

equal  to  Log  C^,  (less  than  Cfa  by  10  ).  The  input  values  of  Temp,,  Log  C^, 
Log  C,  and  Log  C^  are  shown  in  the  inset  of  display  (b)  of  each  figure  along 
with  the  calculated  values  of  Log  n  and  Log  p  ,  the  free  carrier  densities. 

The  calculated  Fermi  energy,  E^,  and  the  screening  factor,  (3,  are  shown 
in  the  inset  of  display  (c).  Also  shown  are  corrected  values  of  the  isolated 
donor  energy,  EQ  ,  (depends  on  0)  and  the  gap  energy,  Ey  (depends  on  T). 


In  the  first  display  of  each  figure,  labeled  (a)  Impurity  Spacing 
Distribution,  the  upper  curve  represents  the  probability  distribution  of  R 
vs  R,  in  units  of  the  isolated  donor  Bohr  radius,  arbitrarily  normalized  to 
a  peak  height  of  unity  for  plotting  convenience.  The  total  area  under  this 
curve  is  proportional  to  the  total  number  of  impurity  pairs  and  is  divided 
into  three  areas,  DD,  DA,  and  AA,  proportional  to  the  number  of  donor-donor, 
donor-acceptor,  and  acceptor-acceptor  pairs,  respectively.  The  areas  between 
these  curves  and  their  general  positions  along  the  R  axis  provide  required 

information  for  weighting  the  energy  levels  vs  R  in  displays  <b)  and  (c). 
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Spacing  between  impurities  is  often  estimated  by  the  rule  d  -  C.  . 

V 

It  is  interesting  to  note  that  this  value  is  too  large  by  almost  85  X,  i.e., 
d  is  larger  than  R  ,  the  most  probable  value  of  R  (peak  of  curve  in  (a)),  by 

IQ 

a  factor  of  1.845,  and  that  the  position  of  d  on  the  curve  is  sufficiently 
far  to  the  right  of  the  peak,  R^,  that  98.5  X  of  the  nearest-neighbor  pairs 
are  more  closely  spaced  than  d.  Only  1.5  X  have  larger  spacing. 

In  compensated  semiconductor  materials  with  concentration  profiles  such 
as  those  shown  in  Fig.  1,  the  donor  and  acceptor  concentrations  are  equal  at 
the  junction.  In  this  case  the  relative  probabilities,  DD,  DA,  and  AA  are 
in  proportion  to  to  ;  te  ;  .  At  points  away  from  the  junction,  DA  numbers 

continue  to  dominate  until  the  concentration  ratio  exceeds  Z  ;  1  (refer  to 
Section  2  and  to  Appendix  B). 

The  part  of  the  probability  distribution  labeled  DD  applies  to  the 
donor  energy  levels,  El  and  EZ,  shown  in  display  (b),  the  first  and  second 
electron  energies  in  a  donor-donor  pair.  The  distribution  labeled  DA  is 
applicable  both  to  E3  and  E4  in  (b)  and  also  to  H3  and  H4  in  (c).  The  AA 
part  applies  to  hole  energy  levels  HI  and  H2  in  (c). 

The  position  of  the  peak  in  the  impurity  separation  distribution  is 
determined  by  the  total  impurity  concentration,  .  As  this 

concentration  is  increased,  from  one  example  to  another,  the  position  of 
the  peak  moves  toward  lower  values  of  R,  thus  changing  the  section  of  an 
energy  level  curve,  in  (b)  or  (c),  that  is  applicable  for  a  given  example. 


Energy  Levels 

Four  levels  are  shown  in  display  <b)  Donor  Energy  Levels.  Of  these, 

El  and  E3  are  the  most  likely  and  nest  most  likely  to  be  occupied  by  a 
localized  electron  since  they  are  single  electron  systems  and  El  is  the 
deeper  of  the  two.  Level  E2  does  not  come  into  existence  unless  level  El 
is  already  occupied  by  an  electron.  Level  E4  comes  into  existence  only  if 
the  DA  pair  is  already  occupied  by  a  hole.  Thus  levels  El  and  E3  are  the 
donor  states  of  primary  importance  in  applications  falling  near  room 
temperature.  Levels  E2  and  E4  become  important  at  low  temperatures  when 
El  and  H3  are  largely  occupied.  A  similar  discussion  applies  to  the  hole 
levels.  H2  and  H4  exist  only  if  HI  and  E3  are  occupied. 

As  the  free  carrier  densities  n  and/or  p  increase,  either  due  to  an 
increase  in  temperature  or  because  of  an  increase  in  the  net  concentration 
C  =  IC^  -  Cfa|  ,  the  value  of  the  screening  factor,  (3,  increases.  This  makes 
the  screening  range,  1/13,  smaller  so  that  neighboring  systems  interact  more 
weakly  and  localized  carriers  are  more  weakly  bound  to  impurity  ions.  In 
other  words,  all  levels  become  shallower,  at  all  values  of  R.  This  trend  is 
clearly  evident  in  the  figures,  especially  at  higher  concentrations. 

In  displays  <a),  (b),  and  <c),  the  units  used  in  the  horizontal  axes 
are  the  same,  the  Bohr  radius  of  an  isolated  donor  electron.  The  hole 
states  have  been  assumed  to  have  larger  Bohr  radii  so  that  the  energy  level 
curves  in  (c)  appear  expanded  in  the  horizontal  direction.  Similarly  the 
unit  of  energy,  effective  Rydberg,  is  the  same  in  both  (b)  and  (c),  but, 
since  the  hole  binding  energy  was  assumed  to  be  smaller,  the  hole  levels 
appear  compressed  in  the  vertical  direction  as  compared  with  the  electron 
levels. 

As  R  increases,  the  levels  E2  and  E4  (also  H2  and  H4)  approach  their 
final  common  value,  E^  (screened  atom  value),  at  much  smaller  values  of  R 
than  in  the  case  of  El  and  E3  (or  HI  and  H3).  At  R  greater  than  two  or  three 
Bohr  radii  there  is  no  further  change  in  the  E2,  E4,  H2,  or  H4  levels.  The 
odd  numbered,  one  electron  or  one  hole  levels  have  not  yet  converged  to  E^ 
even  at  eight  Bohr  radii  in  most  cases.  Since  the  even  numbered  levels  are 


the  most  active  at  low  temperatures  and  the  odd  ones  are  dominant  at  the 
usual  operating  temperatures  of  solar  cells  it  would  seem  that  test  data 
obtained  at  low  temperatures  should  be  projected  with  caution  into  room 
temperature  applications. 

Two  ionized  donor  impurities,  separated  by  a  distance,  R,  form  a 
double  potential  well  with  a  barrier  at  the  midpoint  of  R.  If  the  energy 
of  an  electron,  localized  at  the  well,  is  above  the  potential  barrier,  the 
electron  is  not  restricted  to  one  of  the  ion  sites.  On  the  other  hand,  if 
the  electron  level  is  below  the  barrier  potential,  transport  from  one  site 
to  the  other  requires  a  tunneling  process  involving  time  constants  of 
sufficient  length  to  diminish  contributions  to  conductivity.  The  expression 
for  the  barrier  height  is  given  by  *  -  8  exp(-(3R/2)/R  in  the  usual 
units.  This  function  is  plotted  in  graphic  <b)  as  a  dotted  curve.  It  is 
significant  that,  although  the  donor  energy,  El  (also  E2),  moves  to  lower 
levels  as  R  decreases,  the  barrier  energy,  V^,  moves  to  lower  levels  at  a 
greater  rate,  crossing  El  at  about  six  Bohr  radii  and  moving  substantially 
lower  as  the  ion  separation  becomes  smaller,  This  behavior  can  be  viewed 
as  an  indicator  of  the  onset  of  impurity  conduction.  When  the  bulk  of  the 
donor  pairs  are  closer  than  the  critical  R,  as  shown  by  the  distribution 
DD  in  (a),  the  donor  state  electrons  are  no  longer  strictly  localized.  The 
exact  condition  at  which  this  happens  becomes  a  problem  in  percolation 
theory.  Electrons  with  levels  somewhat  below  the  barrier  level  probably 
contribute  to  the  so-called  hopping  conduction  which  has  a  temperature 
dependence  suggesting  the  tunneling  process.  Note  that  the  curve,  V^, 
shifts  in  position  from  one  example  to  another  because  of  the  dependence 
on  the  screening  factor,  (5. 

The  screening  factor,  13,  also  plays  a  role  in  reducing  the  energy 
gap.  The  energy  needed  to  ionize  a  donor  or  an  acceptor  or  to  form  an 
electron-hole  pair  is  partly  the  result  of  short  range  forces,  such  as 
exchange  interactions  and  ion  core  penetration,  and  partly  due  to  the 
long  range  Coulomb  or  screened  Coulomb  forces.  The  short  range  forces 
are  largely  responsible  for  the  differences  in  the  binding  energies  of 
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different  impurity  ions.  The  additional  contribution  of  the  long  range 
forces  is  essentially  the  same  for  all,  including  the  electron-hole  pair. 

The  calculated  energy  levels  displayed  in  Appendix  A  show  that  the  donor 
and  acceptor  binding  energies  go  to  zero  as  the  value  of  13  increases,  due 
to  the  screening  effect  in  the  long  range  Coulomb  forces.  In  view  of  this 
and  the  remarks  above,  it  follows  that  the  long  range  part  of  the  energy  of 
the  electron-hole  pair  is  also  reduced  by  an  amount  paralleling  the  donor 
energy  reduction  as  (3  increases.  Consequently,  the  energy  gap  will  show 
a  reduction  with  increasing  (3  in  addition  to  changes  for  other  reasons. 

Densities  of  States 

The  density  of  states  distributions  in  displays  (d)  and  (e)  follow 
from  the  probability  distributions  in  (a)  of  the  form,  P(R),  and  the  energy 
levels  in  (b*  and  (c)  of  the  form,  E(R).  A  numerical  routine  is  used  to 
transform  P(R>  and  E(R)  into  the  form  P(E)  which  is  the  required  density  of 
states,  S.  A  separate  transformation  is  performed  for  each  energy  level 
with  the  appropriate  probability  distribution  to  form  displays  <d)  and  (e) 
where  each  distribution  peak  has  the  same  label  as  the  corresponding  energy 
level.  Chart  (d)  contains  the  donor  state  distributions,  along  with  the 
conduction  band  density-of-states  data,  plotted  with  Log  S  as  the  vertical 
axis  and  the  energy,  E,  measured  from  the  intrinsic  conduction  band  edge  in 
effective  donor  Rydbergs,  as  the  horizontal  axis.  Chart  (e>  contains  the 
acceptor  state  distributions  and  the  density-of-states  data  for  the  valence 
band.  The  horizontal  axis  is  in  this  case  referred  to  the  edge  of  the 
intrinsic  valence  band  and  labeled  E  -  E  .  The  units  are  effective  donor 
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Rydbergs.  The  units  of  S  arf  cm  per  effective  Rydberg. 

The  distributions  of  primary  importance,  El,  E3,  HI.  H3,  and  the 
conduction  and  valence  band  states  are  displayed  with  a  solid  line  style. 

The  curves  E2,  E4,  H2,  and  H4  are  shown  with  a  dotted  line.  Additional 
distributions  show,  with  dashed  lines,  filled  states,  i.e.,  those  states 
in  the  bandgap  that  are  occupied  by  localized  electrons  or  holes  and  those 
states  in  the  conduction  or  valence  bands  that  are  occupied  by  carrier 
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electrons  or  holes. 

The  upper  and  lower  portions  of  the  bandgap  are  shown  in  (d)  and  (e) 
with  an  expanded  scale  since  there  is  no  need  for  displaying  the  central 
part.  The  portion  not  shown  varies  between  figures  based  on  different 
temperatures  since  the  value  of  Ev  in  intrinsic  material  depends  upon  the 
temperature.  This  value  is  shown  both  at  the  zero  position  of  (e>,  marked 
EV,  and  in  the  inset  data  in  (d).  In  those  cases  where  the  Fermi  energy 
is  within  the  range  of  E  displayed,  the  position  is  marked  by  a  vertical 
dotted  line  labeled  EF.  In  all  cases  EF  is  given  in  the  inset  data  of  <d) 
referenced  to  the  conduction  band  edge. 

The  shapes  and  positions  of  the  distributions  in  <d)  and  <e)  as  well  as 
their  changes  with  concentration  and  temperature  reveal  quite  clearly  the 
influence  of  neighboring  impurity  interactions  and  of  charge  screening  in 
producing  substantial,  and  in  many  cases  undesirable,  modifications  in  the 
characteristics  of  semiconductor  materials  and  in  diodes  such  as  those  used 
in  solar  energy  conversion. 

One  of  the  most  important  aspects  of  these  effects  is  that  they  are 
inherent  in  the  nature  of  the  materials.  They  are  still  present  in  crystals 
of  the  highest  quality  containing  the  specified  dopants  of  the  greatest 
purity . 

In  the  asymmetric  distribution  of  the  energy  level,  El,  (also  HI) 
of  an  electron  at  a  donor  pair,  (DD>-  ,  there  is  an  impurity  band  tail 
that  extends  toward  the  center  of  the  bandgap  by  values  up  to  four  Rydbergs 
in  the  case  of  donor  pairs.  Clusters  of  three,  four  or  more  donors  were 
not  considered  here  in  order  to  keep  the  calculation  manageable.  In  real 
devices,  such  clusters  can  be  expected,  perhaps  with  lower  probability,  and 
in  those  cases  the  tail  can  be  expected  to  extend  much  farther  into  the  gap. 
Electrons  will  be  localized  in  such  pairs  or  clusters  at  deeper  levels  with 
greater  probability,  thus  causing  them  to  act  as  deep  traps  exactly  as  the 
deep  traps  due  to  undesirable  chemical  impurities. 

In  the  case  of  the  distribution  corresponding  to  the  energy  level,  E3, 
(also  H3)  associated  with  a  donor-acceptor  pair  with  a  localized  electron, 


(DA)-,  the  impurity  band  tail  is  again  present  but  asymmetric  in  the  reverse 
sense  to  El,  with  the  tail  extending  toward  the  conduction  band.  Of  course 
this  reflects  the  fact  that  a  neighboring  acceptor  ion  causes  an  electron 
to  be  bound  more  weakly  to  a  donor  ion.  At  increasing  concentrations,  the 
tail  extends  to  the  conduction  band  and  then  shifts  to  become  a  tail  of  the 
conduction  band  extending  into  the  bandgap.  Near  the  junction,  such  levels 
exceed  the  El  levels  in  frequency  of  occurrence  roughly  by  a  factor  of  two 
and  are  more  likely  to  support  impurity  band  conduction.  Under  these 
conditions,  the  distributions,  E3  and  H3,  are  acting  to  shift  the  effective 
edges  of  both  the  conduction  band  and  the  valence  band  into  the  bandgap, 
thus  narrowing  the  gap  substantially. 

Several  effects  caused  by  changes  in  temperature  are  displayed  in  the 
density  of  states  charts.  In  addition  to  the  inherent  increase  in  bandgap 
with  decreasing  temperature,  present  even  in  intrinsic  materials,  there  is 
a  decrease  in  free  carrier  concentration  which,  in  turn,  causes  a  decrease 
in  (3,  the  screening  constant,  resulting  in  a  general  shift  away  from  the 
band  edges  of  the  distribution  peaks.  In  addition,  as  carriers  become 
localized,  the  one  electron  and  one  hole  levels,  El,  E3,  HI,  and  H3,  become 
filled,  shifting  emphasis  to  the  EZ,  E4,  H2,  and  H4,  second  particle  levels. 
Under  these  conditions,  the  distributions  displayed  with  dotted  lines  are 
the  valid  ones  to  consider  and,  being  much  less  broadened  and  more  nearly 
placed  at  the  unperturbed  dopant  energy  levels,  appear  to  belong  to  a  more 
nearly  ideal  semiconductor  material  with  low  dopant  concentration. 

The  distributions  for  one  particle  states  (solid  line,  odd  numbers)  and 
those  for  second  particles  (dotted  line,  even  numbers)  are  each  plotted  at 
their  maximum  density  values  even  though  the  second  particle  states  exist 
only  to  the  extent  that  the  one  particle  states  are  occupied,  and  the  one 
particle  states  are  diminished  in  number  to  the  extent  that  the  second 
particle  states  are  occupied. 

For  purposes  of  keeping  the  calculation  manageable,  various  device 
circuit  conditions,  e.g.,  open  circuit,  short  circuit,  etc.  were  not 
treated.  The  Fermi  energy,  screening  factor,  and  carrier  densities  were 
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obtained  from  a  thermal  equilibrium  condition  in  bulk  material.  In  the 
space  charge  region  and  in  the  depletion  region  these  parameters  will  be 
modified  slightly  with  the  net  result  that  changes  across  the  junction 
take  place  less  rapidly  and  over  larger  displacements  depending  on  the 
width  of  the  depletion  region,  a  function  of  junction  bias. 
p-n  Fluctuations 

The  calculations  made  in  this  investigation  emphasize  statistically 
random  variations  in  impurity  separations  which  take  place  over  distances 
ranging  from  a  few  A°  up  to  a  few  hundred  A0,  i.e.,  from  the  host  lattice 
constant  up  to  a  few  effective  Bohr  radii.  With  larger  distance  scales  the 
effects  of  random  fluctuations  decrease  in  magnitude,  however,  at  distances 
of  the  order  of  the  width  of  the  depletion  region,  these  effects  may  still 
be  substantial. 
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The  Poisson  relation  P^cv)  =  {(cv>  /k !  >  exp(-cv)  introduced  above, 
which  states  the  probability  of  finding  k  impurities  in  the  volume,  v,  if 
the  average  (mean)  number  is  cv,  has  the  property  of  exhibiting  a  peak  value 
at  k  =  cv  and  also  has  a  standard  deviation,  s  =  (cvp  ,  which  serves  as  a 
measure  of  the  distribution  width.  It  is  useful  here  to  define  a  fractional 
standard  deviation,  f  =  s/cv  =  (cv)  ,  the  ratio  of  the  standard  deviation 
to  the  mean  value.  The  basic  quantities  treated  in  the  relation  are  numbers 
and  their  magnitude  depends  on  size  of  the  volume  under  consideration.  For 


example,  with  c  =  10  cm  ,  if  v  *  10 


then  cv  =  10  and  f  =  1  N  ; 


if  v  =  1C 


,  cv  =  100  and  f  =  10  X  ,  etc. 


The  relations  can  be  put  in  more  familiar  terms  by  calculating  the 
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deviations  in  concentration  in  cm  as  D  =  cf  ,  remembering  that  the  v 

-fc  c 

which  appears  in  f  =  (cv)  is  the  actual  volume  under  consideration  and 


c  is  the  concentration  in  cm 


Then  D  =  (c)^(v)  *  (c/v)* 

c 


For  the  example  shown  in  Figure  1  (a),  near  the  junction  C  and  C 
16  x  b 

are  approximately  10  .  If  the  width  of  the  depletion  region  is  about 
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0.1  microns  then  the  volume,  v,  can  be  taken  to  be  0.001  microns  or 

10  ^  cm^.  Then  D  =  3.16x10^  each  for  C  and  C.  and,  since  deviations 

c  zb 

in  C  and  C.  are  independent,  they  are  large  enough  as  compared  with 


the  net  concentration,  C  =  I -  C^l ,  to  cause  random  fluctuations  in  the 
local  position  of  the  junction  surface  about  the  mean  position  and  produce 
small  closed  regions  or  p-n  anomolies  near  the  junction  in  which  the  p-n 
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type  is  reversed  from  normal.  In  Figure  1  <b),  C  and  C.  are  about  10 

x  b 

near  the  junction.  In  this  case,  for  the  same  v,  the  fractional  s.d.  is 
f  =  10  X  ,  and  =  10*®,  again  larger  than  Cne^  near  the  junction. 

These  types  of  situations  can  be  expected  in  cases  where  the  total 
concentration  is  high  in  the  vicinity  of  the  junction.  The  question  that 
arises  is  what  effects  they  may  have  on  recombination  rates,  efficiency, 
and  power  losses. 


4 .  Summar 


*  - 

|  «  The  model  used  in  calculations  for  this  investigation,  known  as  the  BGN 

■  model,  incorporates  statistically  random  spacing  between  neighboring  dopant 

j,  positions  in  a  GaAs  or  Si  host,  impurity  band  formation  due  to  interactions 

i  of  neighboring  donors  and  acceptors,  screened  Coulomb-like  potential  energy 

|  based  on  the  Thomas-Fermi  approximation  involving  charge  carrier  densities, 

j  effects  due  to  the  higher  concentrations  used  in  compensated  materials,  and 

computing  techniques  designed  to  achieve  accuracy  and  self-consistency  in 
the  process. 

The  results  of  the  calculations,  displayed  in  the  form  of  quantitative 
plots  in  Appendix  A,  show  a  number  of  effects  that  substantially  agree  with 

experimental  knowledge  of  these  semiconductor  systems. 

r 

j  (1)  Donor  and  acceptor  energy  levels  show  broadening  with  increasing 

|  total  dopant  concentration  and  shifting  toward  bandgap  edges  with  increasing 

■  charge  carrier  concentrations. 

(2)  Donor-donor  and  acceptor-acceptor  neighbor  interactions  cause  an 
asymmetrical  impurity  level  broadening  by  introducing  an  impurity  band  tail 

|  reaching  deep  into  the  bandgap. 

(3)  Donor-acceptor  neighbor  interactions  produce  asymmetric  impurity 

'  band  broadening  with  tails  extending  toward  the  bandgap  edges  at  moderate 

l  concentrations.  At  higher  concentrations  and/or  higher  charge  carrier 

l  densities,  the  tail  moves  to  the  bandgap  edge  and  becomes  a  conduction  or 

valence  band  tail  that  is  a  primary  source  of  bandgap  narrowing. 

<4>  Four  sources  of  bandgap  narrowing  are  identified  in  addition  to 
the  temperature  dependent  effect  known  to  be  present  in  the  host  intrinsic 

i 

material.  In  addition  to  the  band  edge  tails  described  in  (3)  above,  the 
electron-hole  pair  creation  requires  less  energy  with  the  screened  potential 
and,  in  the  depletion  region,  electrons  with  energy  levels  lower  than  the 
intrinsic  conduction  band  but  higher  than  the  potential  barrier,  ,  may 


'  '.-VI  ■  S.-W  ^  \Ti  _V  W.  ■«->  V  ^  A  V  V"  *-  -.’  -.  ■.  ••■.  - 


become  mobile,  thus  lowering  the  effective  bandgap  edge.  Holes  behave  in  a 
similar  way  at  the  valence  band  edge. 

(5)  Impurity  band  broadening  effects  are  primarily  due  to  increasing 
values  of  total  concentration,  +  C^,  although  electrical  effects  are 

controlled  by  the  net  concentration,  1C  -  C.  I  ,  in  compensated  materials. 

X  D 

Thus,  broadening  is  greater  in  compensated  materials  than  in  electrically 
similar  uncompensated  materials. 

(6)  The  more  closely  spaced  pairs  or  clusters  of  donors  or  of 
acceptors  have  levels  in  the  bandgap  sufficiently  deep  to  act  as  if  they 
are  a  deep  level  variety  of  chemical  impurity  traps. 

(7)  At  low  temperatures,  with  most  of  the  carriers  localized  and  the 
free  carrier  screening  effect  removed,  the  active  dopant  energy  level 
distributions  appear  sharp  and  unperturbed  as  if  the  impurities  were  dilute 
and  non-interacting,  even  at  relatively  high  concentrations.  At  normal 
operating  (room)  temperatures  or  higher,  different  sets  of  dopant  energy 
levels  become  active  and  screening  effects  are  restored,  resulting  in 
impurity  band  broadening,  bandgap  narrowing,  band  edge  tailing,  etc. 

Low  temperature  data  should  be  projected  into  room  temperature  applications 
with  caution. 

(8)  In  compensated  materials  near  the  location  of  the  P/N  junction, 
statistical  fluctuations  in  concentration  are  sufficiently  large  to  produce 
random  shifts  in  junction  position  about  the  mean  as  well  as  anomolous  cells 
with  inverted  p-n  type  which  may  contribute  to  recombination  rates. 

(9)  The  calculation  reveals  an  anomolous  behavior  at  lowered 
temperatures  resembling  that  of  a  solidstate  phase  transition.  At  a  critical 
temperature  which  varies  with  host  material  and  with  impurity  concentration, 
two  contrasting  stable  conditions  are  possible.  In  one  the  carriers  are 
localized,  screening  is  removed  and  dopant  energy  levels  are  well  separated 
from  the  band  edge.  In  the  other  the  carriers  are  ionized,  screening  is 
restored  and  the  levels  merge  with  the  band  edge.  If  a  device,  such  as  a 
photodiode  in  the  first  state,  received  a  strong  burst  of  radiation,  it 
appears  that  it  would  be  locked  in  the  second  state  until  the  temperature 
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ere  lowered  sufficiently  to  definitely  return  the  system  to  the  first 
itate.  Further  investigation  of  this  type  of  behavior  is  recommended. 
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Fig.  A  101.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El; (DO)-,  E2;(DD-)-,  E3;<DA>-,  E4(DA+>-  vs  R 
(c>  Acceptor  En.  H1;(AA)+,  H2 ; ( AA+ > ♦ ,  H3;(AD)+,  H4(AD->+  vs  R 
(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 


Fig.  A  i02.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b)  Donor  Energy  E1;(DD)-,  E2MDD-)-,  E3;(DA)-,  E4<DA+)-  vs  R 
<c)  Acceptor  En.  Hi>(AA>+,  H2;(AA+>+,  H3>(AD>+,  H4(AD->+  vs  R 
(d>  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  6  H4  vs  E-EV. 
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Fig.  A  103.  (a)  Probability,  AA  >  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  EZi  (DD-)- ,  E3MDA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  H1?<AA>+,  H2 j  ( AA+ ) ♦ ,  H3;(AD>+,  H4<AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  104.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pair*  vs  R 

(b)  Donor  Enargy  El; (DO)-,  E2;(DD->-,  E3; (DA)-,  E4(DA+)-  v*  R 

(c)  Acceptor  En.  H1;(AA>+,  H2)(AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  Statas,  El,  E2,  E3,  &  E4,  vs  E 

(a)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  105.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pair*  v*  R 

(b)  Donor  Enargy  El;<OD)-,  E2;<DD->-,  E3;(DA>-,  E4IDA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;<AA+>+,  H3;(AD)+,  H4(AD-)+  v*  R 

(d)  Dansitias  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  106.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pair*  v*  R 

(b)  Donor  Energy  El;  (DO)-,  E2MDD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;<AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

<e>  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  107.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1MDD)-,  E2;(DD-)~,  E3;(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA)+,  H2;(AA+)+,  H3;(AD)+,  H4(AD-)+  vs  R 
<d>  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  ft  H4  vs  E-EV. 
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111.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

Donor  Eneryy  E1;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
Acceptor  En.  Hi;(AA>+,  H2; < AA+ > +,  H3;(AD>+,  H4(AD->+  vs  R 
Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
Densities  of  Acceptor  States,  HI,  H2,  H3,  6  H4  vs  E-EV. 
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Fig.  A  112.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Ei;(DD)-,  E2;<DD-)-,  E3;(DA)~,  E4IDA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD-)+  vs  R 
Id)  Densities  of  Donor  States,  El,  E2,  E3,  8,  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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113.  <a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

Donor  Energy  E1;<DD>-,  E2;(DD-)-,  E3?(DA)-,  E4(DA+1-  vs  R 
Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  114.  (a)  Probability,  AA  +  DA  ■*•  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2MDD-)-,  E3j(DA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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115.  (a)  Probability,  AA  ♦  DA  +  DD  ,  of  Pair*  vs  R 

Donor  Energy  El ; (DD)-,  E2> (DD- )-,  E3» (DA)-,  E4(DA+>-  v*  R 
Acceptor  En.  Hi;(AA>+,  H2 ;  ( AA+ >  •*■ ,  H3;(AD>+,  H4(AD->+  vs  R 
Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
Densities  of  Acceptor  States,  HI,  H2,  H3,  S  H4  vs  E-EV. 
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Fig.  A  117.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1MDD)-,  E2MDD-1-,  E3;(DA)-f  E4(DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA)+,  H2MAA+1  +  ,  H3;(AD)  +  ,  H4(AD-1+  vs  R 

(d)  Dansitias  of  Donor  Statas,  El,  E2,  E3,  a  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  a  H4  vs  E-EV. 
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Fig.  A  11B.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Elj(DD)-,  E2;(DD-)-,  E3MDA)-,  E4IDA+)-  vs  R 

(c)  Acceptor  En.  Hlj(AA)+,  H2;(AA+>+,  H3>(AD)+,  H4(A0-)+  vs  R 

(d)  Densities  of  Oonor  States,  El,  E2,  E3,  S  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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A  120.  (a)  Probability,  AA  ♦  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;<DD)-,  E2; <DD- )-,  E3>(DA>-,  E4<DA+>-  vs  R 

(c)  Acceptor  En.  H1>(AA)+,  H2;(AA*>+,  H3;(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  124.  (a)  Probability,  AA  ♦  DA  +  DD  ,  of  Pair*  va  R 

<b)  Donor  Energy  E 1 ; ( DD ) - ,  E2MDD-1-,  E3;(DA)-,  E4CDA+1-  vs  R 

(c)  Acceptor  En.  Hi;(AA>+,  H2>(AA+)+,  H3MAD1  +  ,  H4(AD-)+  v»  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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.  A  125.  (b>  Probability,  AA  +  DA  +  DD  ,  of  Pairs 

(b)  Donor  Energy  Elj(DD)-,  E2;(D0->-,  E3;(DA>-,  E4(DA 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3/(AD>+,  H4(AD 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  B  H4  vs 
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Fig.  A  127.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2j(DD-)-,  E3?(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2j(AA+)+,  H3;(AD)+,  H4(AD-)+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  128.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  v*  R 

(b)  Donor  Energy  Elj(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+1-  vs  R 

(c)  Acceptor  En.  H1MAA1  +  ,  H2MAA+1  +  ,  H3;(AD>  +  ,  H41AD-1+  vs  R 
<d>  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  130.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El; (DO)-,  E2;(DD->-,  E3;(DA)-,  E4IDA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA>+,  H2;(AA+>+,  H3;(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  5  H4  vs  E-EV. 
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133.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pair*  vs  R 

Oonor  Energy  El; (00)-,  E2;<00-)-,  E3;(DA>-,  E4(DA*)-  vs 
Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(AD)+,  H4(AD->+  v* 
Densities  of  Donor  States,  El,  E2,  E3,  8>  E4,  vs  E 
Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  138.  (a)  Probability,  AA  +  DA  +  OD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Elj(DD)-,  E2;<DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 


Fig.  A  143.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;<DD)-,  E2MDD-)-,  E3;(LA)~,  E4(DA+>-  vs  R 
<c>  Acceptor  En.  H1;<AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  8>  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 


Fig.  A  14&.  (a)  Probability,  AA  +  DA  *  DO  ,  of  Pairs  vs  R 

<b>  Donor  Energy  E 1 ; <  DD ) - ,  E2MDD-)-,  E3MDA)-,  E4(DA+>-  vs  R 
<c)  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  Ei,  E2,  E3,  6  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  147.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b)  Donor  Energy  Ei;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
<c>  Acceptor  En.  Hi;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  S  H4  vs  E-EV. 
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Fig.  A  148.  <•)  Probability,  AA  +  DA  +  DO  ,  of  Pair*  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2MDD-)-,  E3MDA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  8  H4  vs  E-EV. 
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200.  (a>  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

Donor  Energy  E1?(DD)-,  E2;(DD->-,  E3;(DA)-f  E4(DA+>-  vs  R 
Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD)+,  H4(AD-)+  vs  R 
Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
Densities  of  Acceptor  States,  Hi,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  201.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

(b)  Donor  Enirgy  E1;(DD)-,  E2j (DD- ,  E3;(DA)-,  E4(DA+)-  vs  R 
<c>  Acciptor  En.  H1;<AA>+,  H2MAA+1  +  ,  H3;<AD>+,  H4(AD->+  vs  R 
<d)  Dinsitiis  of  Donor  Statis,  El,  E2,  E3,  &  E4,  vs  E 
<e)  Densitiis  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  204.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

<b)  Donor  Energy  Ei;<DD>-,  E2;(DD-)-,  E3?(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2>(AA+>+,  H3;(AD)+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  20&.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)~,  E2;(DD-)~,  E3MDA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA>  +  ,  H2;(AA+)+,  H3;(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  207.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pair*  vs  R 

(b)  Donor  Energy  El  MOD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  H1;(AA>+,  HZ>(AA+>+,  H3;(AD)+,  H4(AD->+  v*  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  R  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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F i 9 .  A  208.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2MDD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 

( c )  Acceptor  En.  H1,*(AA>  +  ,  H2;(AA+>+,  H3;(AD>  +  ,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States.  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  211.  <a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2?<DD-)-,  E3;(DA)-,  E4<DA+>-  vs  R 
(c>  Acceptor  En.  H1?(AA>+,  H2;(AA+)+,  H3;(AD)+,  H4(AD->+  vs  R 
<d)  Densities  of  Oonor  States,  El,  E2,  E3,  8>  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Z13.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;<DA)-,  E4(DA+)-  vs  R 
Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3; (AD>+,  H4(AD->+  vs  R 
Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  214.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1MDD)-,  EZj(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
<c>  Acceptor  En.  Hi;(AA>+,  H2;(AA+>+,  H3>(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  8  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  8  H4  vs  E-EV. 
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Fig.  A  215.  (a)  Probability,  AA  +  DA  ■»  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
<c>  Acceptor  En.  H1;(AA>+,  H2;(AA+)+,  H3;(AD)+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  216.  (a)  Probability,  AA  +  DA  +  00  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1?(DD)-,  EZj  <DD- )-,  E3;(DA)-,  E4(DA+>-  vs  R 

( c )  Acceptor  En.  H1MAA1  +  ,  H2;(AA+>+,  H3> (AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  Ei,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  217.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;(DA>-,  E4(DA+1-  vs  R 

(c)  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 


.1 

>>.v; 

V  V 


\-y-\V 


s  s 


(a)  Iipuriti)  Spacing  Distribution 


[»I«] 


89 


90 


til  iopuritu  Spicing  Distribution 


R  (lohr) 


R  (lohr) 


(d)  Donor  States  vs  E  (Rgd) 


Cond.  land 


(ei  Receptor  States  vs  E-EV  (Rgd) 


Valence  land 


E-EV  (Rgd) 


Fig.  A  223.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2MDD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  Hi>(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  v-  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  226.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pair*  v*  R 

(b)  Donor  Energy  Elj(DD)-,  E2;(DD-t-,  E3;(DAt-,  E4(DA+t-  v*  R 

( c )  Acceptor  En.  H1;(AA>+,  H2;(AA+t+,  H3;(AD>+,  H4(AD->+  v*  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  K2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  230.  (a)  Probability,  AA  +  DA  ♦  Du  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1MDD1-,  E2MDD-)-,  E3MDA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1/(AA>+,  H2;(AA+)+,  H3;(AD>+,  H4(AD-)+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(el  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  231.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El;  (DO)-,  E2MDD-)-,  E3>  (DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA)  +  ,  H2 ;  ( AA*  >  ■*■ ,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Dansities  of  Donor  Statas,  El,  E2,  E3,  &  E4,  vs  E 

(a)  Densities  of  Acceptor  States,  HI,  H2,  H3,  6  H4  vs  E-EV. 
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Fig.  A  238.  (a)  Probability,  AA  +  DA  DD  ,  of  Pairs  vs  R 

<b)  Donor  Energy  Elj(DD)-,  E2MDD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
<c>  Acceptor  En.  Hi;(AA>+,  H2MAA+1  +  ,  H3;(AD>+,  H4(AD-)+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 
<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  8  H4  vs  E-EV. 


E 


£ 


iii  Iapurity  Spacing  D i sir i but i an 


R  (lehr) 


lb)  Donor  Loyds  (Ryd) 


(c)  Receptor  levels  (Ryd 


TERR. 
LOC  Cb 
IOC  C 
LOC  Cx 
Log  n 
Lay  p 


*  200  K 

*  17.5 
8  17.5 

*  17.801 

*  -18.58S 

*  17.424 


EF-Feroi*  -47.584 
Sern  I  8  0.578 
EB-Sern  8  -0.215 
EV-CRf  *  -48.187 


R  (lehr) 


(0)  Donor  States  vs  E  (Ryd 


Coni,  land 


Y\  El 


(e)  Receptor  States  vs  E-EV  (R^d 


E-EV  (Ryd) 


E  (Ryd) 


»l4s 


HI  Jr  Valence  land 


240.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  v»  R 

Donor  Energy  Elj(DD)-,  EZ»(DD-)-,  E3MDA)-,  E4(DA+>-  vs 
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Fig.  A  241.  (•)  Probability,  AA  +  DA  +  00  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El;  (00)-,  E2;  (DD-)-,  E3MDA)-,  E4(DA+>-  vs  R 
(c>  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(A0)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  244.  (a)  Probability,  AA  ♦  DA  ♦  DD  ,  of  Pairs  vs  R 

<b>  Donor  Energy  Elj(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+>-  v»  R 
(c>  Acceptor  En.  Hi ; ( AA)+ ,  H2 ; < AA+ > ♦ ,  H3MAD1  +  ,  H4(AD->+  va  R 
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(a)  Danaitiaa  of  Acceptor  Stataa.  Hi,  H2,  H3,  ft  H4  va  E-EV. 
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Fig.  A  245.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pair*  vs  R 

(b)  Donor  Energy  Eli (DO)-,  E2;(DD->-,  E3?(DA)-,  E4(0A*)-  v*  R 
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Fig.  A  247.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pair*  va  R 

<b)  Donor  Energy  El; (DO)-,  E2;<00-)-,  E3;(0A)-,  E4(DA+)-  v«  R 
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Fig.  A  248.  (a)  Probability,  AA  ♦  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Oonor  Energy  El  MOD)-,  E2MDD-)-,  E3;(DA>-,  E4(DA+>-  vs  R 
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(d)  Oensities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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302.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 
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Fig.  A  305.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b>  Donor  Energy  E1;(DD>-,  E2;(DD-)-,  E3;(DA)-,  E4<DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  6  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  306.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b>  Donor  Energy  Ei;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+>-  vs  R 
(c)  Acceptor  En,  Hij(AA)+,  H2MAA+1  +  ,  H3; (AD)+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  308.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pair*  v*  R 

(b)  Donor  Energy  El; (DO)-,  E2;<DD-)-f  E3;(DA)-,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  Hi;CAA)+,  H2;(AA+>+,  H3;(AD)+,  H4(AD->+  v*  R 
(<j)  Densities  of  Donor  States,  El,  E2,  E3,  8  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  6  H4  vs  E-EV. 


Fig.  A  311.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Elj(DD)-,  E2;(DD->-,  E3;(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  313.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b>  Donor  Energy  E1;(DD)~,  E2MDD-)-,  E3?(DA)-r  E4<DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  31&.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Elj(DD)-,  E2MDD-)-,  E3MDA)-,  E4(DA+)-  vs  R 
<c)  Acceptor  En.  H1>(AA1+,  H2MAA+1  +  ,  H3MAD1  +  ,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  317.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pairs  vs  R 

<b)  Donor  Energy  El;(DD)-f  E2;(DD->-»  E3;(DA)-,  E4(DA+1-  vs  R 
(c>  Acceptor  En.  H1;(AA)+,  H2MAA+1  +  ,  H3;(AD>+,  H4(AD->*  vs  R 
(d>  Dansities  of  Donor  Statas,  El,  E2,  E3,  B.  E4,  vs  E 
(a)  Dansities  of  Acceptor  States,  HI,  H2,  H3,  8>  H4  vs  E-EV. 
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Fig.  A  324.  (a)  Probability,  AA  +  DA  4  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;<DD-)-,  E3;(DA>-,  E4(DA+t-  vs  R 

(c)  Acceptor  En.  H1»(AA)+,  H2?(AA+>+,  H3j<AD>+,  H4(AD->+  vs  R 
(dt  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  327.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1»(DD)-,  E2;(DD-)-,  E3; (DA)-,  E4<DA+>-  vs  R 

(c)  Acceptor  En.  H1;(AA)*,  H2;(AA+)+,  H3;(AD)+,  H4(AD->+  vs  R 
Id)  Densities  of  Oonor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States.  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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ig.  A  328.  (a)  Probability,  AA  +  DA  ♦  OD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El ?  (DO)-,  E2MDD-1-,  E3;(DA)-,  E4(DA*>-  vs  R 

(c)  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Oonor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  6  H4  vs  E-EV. 
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Fig.  A  331.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b>  Donor  Energy  E1;(DD)-,  E2;<DD->-,  E3;(DA)-,  E4IDA+)-  vs  R 

(c)  Acceptor  En.  H1>(AA)  +  ,  H2MAA+>  +  ,  H3;(AD)  +  ,  H4IAD-)-*-  vs  R 

(d)  Densities  of  Oonor  States,  El,  E2,  E3,  &  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  337.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Oonor  Energy  El; (DO)-,  EZ;(DD->-,  E3; (DA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  Hl;(AA)-r,  H2;(AA+>+,  H3; (AO>+,  H4(AD->+  vs  R 
<d)  Oensities  of  Oonor  States,  El,  E2,  E3,  R  E4,  vs  E 

(e)  Oensities  of  Acceptor  States,  Hi,  H2,  H3,  R  H4  vs  E-EV. 
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Fig.  A  338.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b )  Donor  Energy  E1;(DD>-,  EZ;<DD->-,  E3;(DA)-,  E4(DA+1-  vs  R 
(cl  Acceptor  En.  Hi;(AA>+,  H2MAA+1  +  ,  H3;CAD>  +  ,  H4(AD-)+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  HZ,  H3,  &  H4  vs  E-EV. 
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Fig.  A  340.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD>-,  E2MDD-)-,  E3?(DA)-,  E4CDA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2j(AA+)+,  H3;(AD)+,  H4<AD-)+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  341.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El;(DD)-f  E2;(DD-)-,  E3;<DA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  H1;(AA>+,  H2;<AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  R  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  R  H4  vs  E-EV. 
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343.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

Donor  Energy  Ei;(DD)~,  E2;<DD-)-,  E3MDA)-,  E4(DA+)-  vs  R 
Acceptor  En.  Hi;(AA>+,  H2;<AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
Densities  of  Donor  States,  Ei,  E2,  E3,  &  E4,  vs  E 
Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  345.  <a)  Probability,  AA  +  DA  +  Du  ,  of  Pairs  vs  R 

(b)  Oonor  Energy  Elj(DD)-,  EZ;(OD->-,  E3;(DA)-,  E4(DA-*,>-  vs  R 
<c)  Acceptor  En.  H1>(AA)  +  ,  H2MAA+)  +  ,  H3»(AD)+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  R  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  346.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pair*  vs  R 

<b)  Donor  Energy  E1;(DD)-,  E2; (DD- )-,  E3; (DA)-,  E4(DA*>-  vs  R 
<c)  Acceptor  En.  Hi;(AA>+,  H2;(AA+)+,  H3;(AD>+,  H4(AD-)+  v»  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  S  H4  vs  E-EV. 
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Fig.  A  401.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;<DD-)-,  E3;(DA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  Hi;(AA)+,  H2;(AA+>+,  H3;(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV, 
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Fig.  A  402.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD>-,  E2MDD-)-,  E3;(DA)-r  E4<DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA>+,  H2;<AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  435.  (•)  Probability,  AA  ♦  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El; (DO)-,  E2;(DD-)-,  E3;<DA)-,  E4<DA+)-  vs  R 
<c)  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(AD)+,  H4(AD-)+  v»  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  S  H4  vs  E-EV. 
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<b)  Donor  Energy  E1;(DD)-,  EZj(DD-)-,  E3j(DA>-,  E4(DA+)-  vs  R 
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(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  R  H4  vs  E-EV. 
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i  A  442.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;<DD)-,  E2M0D--1-,  E3;(DA)-,  E4<DA+>-  vs  R 
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(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  443.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El;  (DO)-,  E2MDD-)-,  E3;<DA)-,  E4(DA+)~  vs  R 

(c)  Acceptor  En.  Hi;(AA>+,  H2;<AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  444.  <a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Elj(DD)-,  EZj(DD-)-,  E3;(DA)-,  E4<DA+>-  vs  R 

(c)  Acceptor  En.  H1?(AA)+,  H2j(AA+)+,  H3;(A0>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(#)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  447.  (a)  Probability,  AA  ♦  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Ei»(DD)-,  E2;<DD-)-,  E3»(DA)-,  E4(DA+>-  vs  R 
(c>  Acceptor  En.  H1;(AA>+,  H2 / < AA+ ) ♦ ,  H3;(AD>+,  H4(AD-)+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
<e>  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  500.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;(0A>-,  E4IDA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA>+,  H2;(AA+)+,  H3;(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  501.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;<AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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,  A  503.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Oonor  Energy  E1;(DD)-,  E2MDD-)-,  E3;(DA)-r  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  505.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;<DD>-,  E2;<DD->-,  E3;(DA)-,  E4<DA+)-  vs  R 
<c>  Acceptor  En.  HiHAA)*,  H2;<AA+)+,  H3>(AD>+,  H4(AD-),r  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  ft  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  8<  H4  vs  E-EV. 
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Fig.  A  50S.  (a)  Probability,  AA  ♦  DA  ♦  DO  ,  of  Pair*  vs  R 

(b)  Donor  Energy  E1MDD)-,  E2j<DD-)-,  E3;(OA)-,  E4(DA+)-  vs 
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<d)  Densities  of  Donor  States,  El,  E2,  E3,  S  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 


Fig.  A  507.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-r  E3;(DA)-r  E4(DA+)-  vs  R 

(c)  Acceptor  En.  Hi;(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  508.  (a)  Probability,  AA  +  DA  ♦  DO  ,  of  Pair*  v*  R 

(b)  Donor  Energy  El; (DO)-,  E2;(DD-)~,  E3;(DA)-,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  H1;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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,  A  511.  <a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)~,  E2;(DD-)~,  E3;(DA)-,  E4(DA+)-  vs  R 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->-*-  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  512.  (3)  P>~obabi  1  ity,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(t)  Donor  Energy  Elj(DD)-,  E2;(DD->”»  E3;(DA>-,  E4<DA+)-  vs  R 
(c>  Acceptor  En.  H1;(AA)+,  H2;(AA+)+,  H3;(AD>+,  H4(AD->*  vs  R 
(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  514.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b>  Donor  Energy  E1;(DD>-,  E2;(DD->-,  E3/(DA>-,  E4(DA+>-  vs  R 
<c>  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
(d>  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States.  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  515.  (a)  Probability,  AA  ♦  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El; (DO)-,  E2;(DD-)-,  E3j(DA)-f  E4(DA*)-  vs  R 

(c)  Acceptor  En.  Hlj(AA)+,  H2 ;  < AA+  > ♦ ,  H3;(AD)+,  H4(AD->*  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  S  E4,  vs  E 

<e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  6  H4  vs  E-EV. 
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Fig.  A  516.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  El; (DO)-,  EZj(DD-)-,  E3;(DA>-,  E4(DA*)-  vs  R 
(c>  Acceptor  En.  H1;(AA)+,  HZ;<AA*>+,  H3;(AD)+,  H4(A0->*  vs  R 
(d>  Densities  of  Donor  States,  El,  EZ,  E3,  6  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  HZ,  H3,  &  H4  vs  E-EV. 
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Fig.  A  518.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b)  Donor  Energy  E1?(DD)-,  E2;(DD-)-,  E3;(DA)-,  E41DA+)-  vs  R 
(c)  Acceptor  En.  Hi»(AA>+,  H2MAA+)  +  ,  H3;(AD>+,  H4(AD->+  vs  R 
(d>  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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.  A  522.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b>  Donor  Energy  E1;(DD)-,  E2; (DD- )-,  E3?(DA)-,  E4(DA+>-  vs  R 

(c)  Acceptor  En.  H1;(AA>  +  ,  H2;<AA+>+,  H3;(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  6  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV, 
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(a)  Iipurity  Spacing  Distribution 


Fig.  A  52£.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

<b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  Hi;(AA)+,  H2;(AA+)+,  H3?(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  HZ,  H3,  6  H4  vs  E-EV. 
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Fig.  A  528.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b>  Donor  Energy  El; ( DD ) - ,  E2; <DD- E3; (DA)-,  E4(DA+>-  vs 
(c>  Acceptor  En.  H1»(AA>+,  H2;(AA+>+,  H3;(AD)+,  H4IAD-1+  vs 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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,  A  530.  <a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD>-,  E2MDD-)-,  E3;(DA)~,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  H1;(AA>+,  H2 ; < AA+ ) ♦ ,  H3?(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 


E-EV  ((yd) 


Fig.  A  531.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b>  Oonor  Energy  E1;(DD)-,  E2;(D0-)-,  E3;(DA)-,  E4<DA+)-  vs  R 
(c>  Acceptor  En.  H1;(AA>+,  H2;(AA+>+,  H3;(AD)+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  536.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pair*  vs  R 

<b)  Donor  Energy  E1;(DD)-,  E2;(DD->-,  E3;(DA)-,  E4(DA+)-  vs  R 
(c>  Acceptor  En.  Hi>(AA>+,  H2;(AA+>+,  H3MAD)*,  H4(AD-)+  v*  R 
<d)  Densities  of  Donor  States,  Ei,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  Hi,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  540.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  Ei;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
(ci  Acceptor  En.  Hi;(AA>+,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
<d)  Densities  of  Donor  States,  El,  E2,  E3,  R  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  8  H4  vs  E-EV. 
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(a!  Iipuritg  Spacing  Distribution 


Fig.  A  542.  (a)  Probability,  AA  +  DA  +  DO  ,  of  Pairs  vs  R 

<b)  Donor  Energy  E1;(DD)-,  E2;<DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
(c)  Acceptor  En.  H1;(AA)+,  H2?(AA+)+,  H3;(AD)+,  H4(AD-)+  vs  R 
(d!  Densities  of  Oonor  States,  El,  E2,  E3,  &  E4,  vs  E 
(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  8.  H4  vs  E-EV. 
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Fig.  A  544.  (a)  Probability,  AA  +  DA  ♦  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  EljtDD)-,  EZj(DD-)-,  E3;(DA)-,  E4(DA+>-  vs  R 
(c>  Acceptor  En.  HiMAA)*,  H2j(AA+>+,  H3?(AD>+,  H4(AD-)+  vs  R 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  6  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  ft  H4  vs  E-EV. 
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Fig.  A  547.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2?(DD-)-,  E3;(DA>-,  E4<DA+>-  vs 

(c)  Acceptor  En.  H1;(AA)+,  H2;(AA+)+,  H3;(AD)+,  H4(AD-)+  vs 

(d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E 

(e)  Densities  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Fig.  A  548.  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R 

(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3;(DA)-,  E4(DA+)-  vs  R 
( c  >  Acceptor  En.  H1>(AA)  +  ,  H2;(AA+>+,  H3;(AD>+,  H4(AD->+  vs  R 
(d)  Dansitias  of  Donor  Statas,  El,  E2,  E3,  &  E4,  vs  E 
(a)  Dansitias  of  Acceptor  States,  HI,  H2,  H3,  &  H4  vs  E-EV. 
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Donor-Donor  Pair  Energy  Calculations 


The  procedure  makes  use  of  a  variation  type  calculation  similar  to  that 
used  by  Slater  in  calculating  the  ground-state  energies  of  H2  and  H2+. 

The  impurity  centers  are  located  at  points  a  and  b  separated  by  a  distance  R. 
Distances  from  these  centers  to  electrons  at  points  1  and  2  are  denoted  by 
rla»  rlb>  r2a»  an<*  r2b*  The  inter-electron  distance  is  r^2.  Two 
angles  4>^,  and  locate  the  azimutal  positions. 

3  1/2 

Four  normalized  Is  wave  functions  of  the  form  a(l)  *  (a  /ir)  exp(-ar,  ) 

la 

are  used  to  form  the  system  wave  function  along  with  appropriate  spin 
functions.  For  the  H2  type  system,  the  ground  state  function  has  the  form: 

*UA  =  [a(l)b(2)  +  a(2)b(l)]UA  (1) 

where  is  an  antisymmetric  spin  function.  The  spin  function  does  not 
enter  into  the  energy  integral  but  determines  the  symmetry  of  the  orbital 
function  for  the  ground  state  (the  Hamiltonian  is  spin  independent).  In  the 
H2+  case  the  -orbital  function  has  the  form: 

'P  -  a(l)  +  b(l)  (2) 

The  Hamiltonian  for  the  H2  type  system  (in  units  of  Bohr  radii  and 
Rydbergs)  has  the  form: 

2  2 
H  *  -V.  - 


-  V. 


-  V..  -  V. 


-  V_.  +  V.  - 


(3) 


where  the  potential  energy  terms  are  of  the  form  Vla  -  2exp(-£ria) /r^a 

> 

etc.,  and  6  is  the  inverse  screening  range.  For  the  H2+  system,  terms 
with  subscript  2  are  omitted. 

The  integrals  to  be  calculated  come  from  the  expression 

E  ■  Jty1'  H  ip  dv/  jV*dv  (4) 

for  each  of  the  two  systems.  The  parameter  a  in  the  basis  Is  function  is  the 
variation  parameter.  It,  of  course,  has  the  value  unity  for  an  isolated 
hydrogen-like  atom,  but,  in  the  case  of  a  molecule-like  system  it  is  adjusted 
by  minimizing  the  energy  with  respect  to  a  variation  of  a.  The  above 
integrals  can  be  expressed  in  terms  of  the  following  integrals  over  the  basis 


functions:  (w  -  aR) 

y*a(l)  (-V12)a(l)dv1  -  a2  (5) 

J”a(l)b(l)dv^  *  S  =  e  W<1  +  w  +  w2/3)  (6) 

/a(l)(-V12)b(l)dv1  =  -a2 (K  +  S)  =  a2e"w(l  +  w  -  w2/3)  (7) 

/a2(l)(-Vla)dv1  -  -2a/ (1  +  S/2a)2  (8) 

/a2(l)(-Vlb)dv1  -  ajg  (9) 


/aUmiX-V^JdVj,  -  aKg 


(10) 


a2(l)b2(2)V 


12dvldv2 


aj 


I 


(11) 


/a(l)b(l)a(2)b(2)V12dv1dv2  -  ctKg  (12) 

Other  Integrals  encountered  In  the  evaluation  of  Equation  (4)  can  be  reduced 
to  those  above  by  interchanging  indices.  Expressions  for  the  first  four 
relatively  simple  integrals  are  included  above.  Expressions  for  Equations 
(9),  (10),  and  (11)  are  long  and  cumbersome  and  will  be  treated  in  another 
section.  Equation  (12)  cannot  be  integrated  in  closed  analytical  form  and 
requires  numerical  evaluation. 

The  first  three  of  the  above  Integrals  do  not  depend  explicitly  on  g 
in  contrast  to  the  remaining  five.  The  last  four  are  functions  of  a,  3,  and 
R,  however,  if  the  substitutions  g  -  3/2CX  and  w  -  aR  are  made,  then  Jg, 

Ks,  j'3  and  k'8  are  functions  of  g  and  w  only.  In  the  process  of 
minimizing  the  energy  by  varying  a,  the  final  value  of  a  depends  on  3  so  that 
the  first  three  integrals  depend  implicitly  on  0.  The  last  five  depend  both 
explicitly  and  implicitly  on  3.  The  numerical  evaluation  of  x's,  to  be 
treated  in  detail  later,  covered  an  array  of  1025  values,  with  0  <  g  _<  2.4  in 
steps  of  0.1,*  and  0  w  <  8.0  with  steps  of  0.2.  A  series  of  the  form 
G(g)e~2w(l  +  A^w  +  A2w2  +  A3W3  +  A^w4)  with 

G(g)  -  (f2/2  +  f3  +  f4)/2,  was  fitted  to  this  array  by  least  square  ad¬ 
justment  of  the  A's  for  each  g,  good  to  four  decimal  places  (f  *  1/(1  +  g) ) • 
In  terms  of  the  above  expressions,  the  energy  becomes,  with 
f  -  (1  +  3/2a)”l, 

Ex(a)  «  (1  +  S)-1(a2(l  -  S  -  K)  +  a(-2f2  +  Jg  +  2Kg))  (13) 


for  the  H2+  -  type  system,  and 


for  the  two  electron  il2  -  type  system. 

These  expressions  must  be  minimized  independently  with  respect  to  a 
variation  of  for  each  value  of  3  and  R  selected.  This  procedure  for 
calculating  molecular  binding  energies  has  been  shown  to  give  results  within 
2  or  3  percent  of  correct  values  over  the  mid-range  of  R  and  even  better 
accuracy  at  large  and  small  R. 

This  calculation  leads  to  two  binding  energies  of  interest;  first 
the  energy  of  one  electron  bound  to  a  pair  of  impurities  at  separation  R  and 
concentration  c,  second  E2  ■  S^2  “  E]_  the  energy  for  binding  a  second 
electron  to  the  same  system. 


Jm J  ». 


Density  of  States  for  and 

Because  of  the  variation  of  the  binding  energies,  E^  and  E2,  with 
impurity  spacing,  R,  and  with  screening,  g  ,  the  donor  level  (or  acceptor 
level)  cannot  be  considered  as  a  constant,  independent  of  impurity 
concentration.  The  resulting  band  of  levels  can  be  described  by  two 
distributions,  one  for  E^  and  one  for  E2,  which  can  be  calculated  from 
E^  and  E2  (as  functions  of  R  and  c),  and  P(R,  c),  where  c  here  is  the 
total  concentration,  i.e.,  c  =  +  C^. 

The  value  of  P(R,  c')  is  given  by 

P(R,  c')dR  =  (4irc')R^dR  exp(-4irc'R3/3)  (15) 

where  R  is  in  effective  Bohr  radii  and  c'  is  concentration  in  corresponding 
units.  If  c  is  the  concentration  in  (cm)~^,  then,  the  frequency  of  pairs 
with  separation  between  R  and  R  +  dR  per  cm-*  is 

N(R)dR  =  (c/2)P(R,  c')dR  =  (c/2)PR  (16) 


From  the  values  of  E^  (say)  at  R  and  at  R  +  dR  (where  dR  is  small  but 
finite)  the  corresponding  values  of  E^  and  E^  +  dE]^  are  determined. 
Consequently,  the  frequency  in  Equation  (16)  is  identified  with  the  frequency 
falling  between  E^  and  E^  +  dE^,  thus 

S1(E1)  =  (c/2)P(R,  c')dR/dE1  (17) 
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All  of  Che  quantities  on  the  right  side  are  known  at  least  in  numerical  form 
so  that  a  plot  of  S^(E^)  and  similarly  82(^2)  can  be  produced  with  an 
appropriate  computer  program.  Care  must  be  exercised  In  implementing  this 
program  because  of  the  behavior  of  dE^  and  dE2  which  may  occasionally 
become  zero. 

The  calculation  produces  two  overlapping  band-like  distributions  of 
donor  levels  (or  acceptor  levels)  which,  upon  increasing  concentration, 
broaden  and  form  a  "tail"  which  extends  deep  into  the  band  gap.  There  is  an 
accompanying  decrease  in  the  energy  separation  between  the  conduction  band 
and  the  peak  of  the  impurity  band  distribution. 


.  '  .  V .  >■/ 


,v.  •  .V.V.VlW.V V. V  V  s' 


Recalculation  of  Fermi-level  and  Screening  Factor 
The  starting  values  for  the  screening  parameter,  g,  and  the  Fermi-level, 
Ep  »  kTriy  are  obtained  from  the  charge-balance  equation 


a  *  N+  +  p  (donor  case) 


and  the  relation 


B2  =  (q2/ekT)dn/dru 


based  on  the  Thomas-Ferai  approximation.  The  conduction  electron 
concentration,  n,  has  the  form 


n  -  21f1/2NcF1/2(nF) 


with  Nc  the  effective  density  of  states  in  the  conduction  band  and 
F^/2(Hf)  the  Fermi  integral 


00  1  /  2  —1 

Fi/2(r¥}  =  /  n  (1  +  ex^(n  ‘  V>  dn 


Equation  (20)  can  be  put  in  the  form 


n  -  Nc  exp(nF)SF 


m^*>~**  *-V  jc-V  ,«.A.  j>.A  mjf .  a  -  a*.  —  **-  *  --  i*-  V—  *"»-•  -  »*-  ^  *'  * 


-r  ^  ^  v  »■ ru*  _■  w*  -■  v  ^  ^  ">  -j*  >  wj<  ”-¥  -j,  J*  s  r  >  .»  ,  .' '>  V  7  fTs  7  /  / 


where  Sp  has  the  form 


CO 

Sp  =  (tt1^2/2)  J"  n1^2  df|/ (exp  Tip  +  expn), 


The  factor  Sp  has  the  effect  of  correcting  the  Boltzman  statistics,  which 
applies  in  the  expression  n  ■  Ncexp(r)p)  for  the  non-degenerate  case,  to 
Fermi  statistics  for  the  general  case.  A  similar  expression  can  be  developed 
for  p.  In  the  initial  calculation  the  expression  for  N+  is 


N  -  Nd(1  -  F (E) ) 

N+  -  Nd(1  +  g  exp((EF  -  E^/kT))"1 


where  Nq  *  c  is  the  number  of  donor  impurities  and  Ep  is  the  donor  energy 
level,  assumed  to  be  constant.  The  spin  degeneracy,  g,  is  2  for  this  case. 
Equation  24  is  valid  for  low  concentrations  where  the  impurities  are  isolated 
from  one  another  and  the  donor  (acceptor)  energy  level  is  unperturbed. 

At  the  medium  and  high  impurity  concentrations,  the  ground  state 
energies,  E^(R)  and  E2(R) ,  as  calculated  above,  show  substantial 
variation  from  Eg  (as  much  as  4Eg)  and  vary  from  site  to  site. 

Consequently,  the  expression  for  If4-  takes  the  form  of  a  sum  over 
contributions  with  different  energies.  Since  the  energies  are  known  as 
functions  of  R  as  well  as  the  frequencies  of  occurance,  N(R) ,  it  is 
convenient  to  form  the  sum  over  incremental  values  of  R  sudh  that 
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where  dR  is  small  but  finite,  is  the  number  of  impurity  pairs  with  separation 
between  R  and  R  +  dR.  Among  these  Ng  pairs  are  02  neutral,  n^  singly 
ionized  (charge  -  +1),  and  n0  doubly  ionized  (charge  •  +2)  for  the  donor 
case.  From  statistical  mechanics  the  Fermi  statistics  relations  are 


nQ/“i  ■  gx  exp((E1  -  EjJ/kT)  -  P1(R) 

ni^n2  “  82exp^E2  ~  *  P2(r) 

\  "  Qo  +  nl  +  n2 

which,  along  with  the  charge  relation 


(26) 

(27) 

(28) 


0*n2  + 


l*n.  +  2*n 
1  o 


(29) 


and  Equation  (16),  lead  to 


N+  -  c  £r  Pr(1/2  +  pj)/(l  +  Px  +  1/P2) 


(30) 


The  factors  g^  and  g2  are  spin  degenerate  values  and  are  1/2,  and  2 
respectively  for  the  present  case.  Equation  (30)  is  evaluated  numerically 
from  .'JR  »  N(R)dR,  calculated  earlier,  and  tabular  values  of  E^(R)  and 
E2(R),  also  from  the  previous  calculation. 
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This  new  value  for  N*  is  inserted  in  Equation  (18)  along  with  the 
expression  for  n  from  Equation  (22),  leading  to  a  value  for  exp(Ep/kT)  by 
means  of  successive  approximations  of  Equation  (18)  by  an  iterative 
procedure . 

The  new  values  of  Ep  and  n  are  then  used  to  obtain,  from  Equation 
(19),  a  new  value  for  the  screening  parameter,  8.  At  this  point,  if  8  changes 
noticably  from  the  previous  value,  the  values  of  Ep(R)  and  E2(R)  must  be 
redone  and  this  overall  procedure  iterated  until  a  self-consistent  set  of 
quantities  results, i.e. ,  Ep,  n,  8,  Ep(R),  E2(R).  This  must  be  done,  of 
course,  for  each  selected  value  of  concentration  and  temperature. 

Distributions  of  donor  (acceptor)  levels  that  are  occupied  by  localized 
electrons  (holes)  are  calculated  after  the  above  conditions  are  met.  These 
come  from  procedures  similar  to  those  used  for  Si(Ep)  and  S2(E2), 

Equation  (17),  except  that  the  quantities  p^  and  P2  (Equations  26-27)  are 
incorporated  in  the  numerical  calculation. 

In  compensated  materials,  the  expression,  PR,  appearing  in  Equation  30, 
above,  is  multiplied  by  the  factors,  PE  in  the  case  of  donor-donor  pairs, 
and  PD  in  the  case  of  donor-acceptor  pairs,  where  P  =  (l-q)~  ,  P  =  2a(l-q), 
and  q  -  C^/(,C^+C^).  This  serves  to  adjust  the  total  impurity  concentration, 
c  ~  ^2 ’  appropriate  value  for  the  type  of  impurity  pair 


under  consideration. 


In  this  procedure,  two  systems  are  considered,  first  a  donor-acceptor 
impurity  pair  with  a  single  localized  carrier,  either  an  electron  or  a 
hole,  and  second  the  pair  with  both  localized  carriers  present.  The  first 
system  has  a  groundstate  energy  called  and  the  second  From  these 

the  quantity  ~  is  obtained  representing  the  energy  needed 

to  remove  one  of  the  two  from  the  system.  In  the  first  case  it  is  con¬ 
venient  to  fix  attention  on  an  electron  bound  to  a  donor  impurity  with 
an  acceptor  ion  located  at  a  distance  R.  As  before,  donor  ion  to 
electron  distance  is  r^  and  the  acceptor  ion  to  electron  is  In 

this  case  the  center  of  the  isotropic  electron  cloud  (origin  of  the  wave 
function)  is  assumed  to  be  offset  from  the  donor  ion  by  a  small  distance  L 
along  the  direction  of  the  line  of  centers  of  the  donor-acceptor  pair, 
opposite  from  the  acceptor  ion  being  positive.  This  step  serves  to  allow 
for  polarization  of  the  neutral  donor  impurity  by  the  acceptor  ion,  thus 
improving  the  approximation  to  the  system  energy,  but  complicating  the 
computation  by  introducing  a  second  variation  parameter.  The  energy  is 
minimized  with  respect  to  both  the  offset,  L,  and  the  inverse  Bohr  radius, 
a,  thus  adjusting  both  the  polarization  and  the  size  of  the  electron 
cloud.  The  distance  from  the  offset  origin  to  the  electron  is  r^,  so 
that,  in  the  Hamiltonian,  (Equation  3  with  terms  omitted  having  subscript 
2)  the  first  (kinetic  energy)  term  involves  differentiation  with  respect 

to  r, ,  the  second,  -V.  ,  remains  the  same  and  the  third,  +V..  ,  has  the 
1  la  id 

sign  reversed  since  the  acceptor  ion  is  negative.  With  V  =2  exp(-6r  )/ 

13  13 


rmm'-y v ■>  v  ■  ■  ■.•■rv.  '.1'  ? ’> b.- v> *>  v 


H  =  V 


1 


V.  +  V,v 
la  lb 


(31) 


The  wave  function  assumed  for  this  case  is  the  (Is)  basis  function  centered 
at  the  offset  origin. 


3  % 

=  (a  /it)  expt-ar^ 


(32) 


The  resulting  expression  for  the  energy  before  it  is  minimized  can  be  stated 
in  terms  of  the  integrals  (Equations  5-12)  developed  earlier,  except  that 
they  become  functions  of  L  and  (L+R)  rather  than  R  as  in  the  earlier  case. 


E  =  a  +  ctJ  (L)  -  aJ  (L+R) 
1  s  s 


(33) 


In  case  both  the  electron  and  the  hole  are  present,  procedures  similar 
to  those  above  are  used.  The  wave  functions  of  both  particles  are  assumed 
to  have  their  origins  offset  by  equal  magnitudes  and  in  opposite  directions 
along  the  line  of  centers  of  the  impurity  ions.  The  modified  Hamiltonian 
for  this  system  is 


H  =  -V 


1 


-  Vla  +  Vlb  +  V2a  '  V2b  "  V12 


(34) 


and  the  assumed  wave  function  consists  of  a  product  of  the  ls(r^)  and 
the  ls(r2)  basis  functions  centered  at  the  offset  origins  near  the  donor 


ion  and  the  acceptor  ion,  respectively. 


ip  =  (a  /tt )  exp(-ar^)  exp(-ar2) 


(35) 


Again,  the  expression  for  the  energy  before  minimizing  can  be  written 
in  terms  of  the  integrals  (Equations  5-12)  as  functions  of  L,  (L+R),  and 
( 2L+R) , 
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E12  =  2a  +  a(2Jg(L)  -  2Jg(L+R)  -  Jg'(2L+R))  . 


(3 


The  minimized  values  of  the  expressions  in  Equations  33  and  36  are  the 
energies  of  interest,  and  D^,  from  which  D2  =  D12  -  is  also  found 
In  Appendix  A,  these  energies,  as  well  as  those  for  donor-donor  pairs, 
are  plotted  as  a  function  of  R  for  an  array  of  temperatures  and  concen¬ 


trations. 
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APPENDIX  C 


/a2(l)(-V  )dV  =  aJ  where 
J  lb  1  s 


aJg  =  2ay (e  2w  -  y(e  -  e-2w)/w) 


and  y  =  (1  -  g2)  2 
Let  Y  =  (2a  -  B)/R 
Yx  -  (1  -  e"Y)/Y 
Y2  =  (1  -  Y  -  e~Y)/Y2 
Y3  =  (1  -  Y  +  .5Y2  -  e-Y)/Y3 
Y4  =  (1  -  Y  +  .5Y2  -  Y3/6  -  e~Y)/Y4 


(C.5) 


Then 


Y..  =  Y-Y.  +  1/6 
3  A 

Y2  =  Y-Y3  -  1/2 
Y1  =  Y*Y2  +  1 


By  using  the  above  definitions,  the  expression  reduces  to 
-$R, 


aJs  =  -2afe  ^(2w(2fY2  +  Y  )  +  f) 
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/a(l)(-V1Z)a(l)dV1  =  a  2 

(C.l) 

/a(l)b(l)dV1  =  S  •  e_w(l  +  w  +  w2/3) 

(C.2) 

/a(l)(-V12)b(l)dV1  =  -a2(K+S)  =  a2e~W(l+w-w2/3) 

(C.3) 

i.  A  a  „-d 

/ a2(l)(-Vla)dV1  =  -2a (1  +  B/2a)~2  =  -2af2 

(C.4) 

• _ i 

.•V 

where  f  *  (1  +  g)  1,  g  =  B/2a 

vV  /'J 
-.'••■.•A 

*  * . ' .« .i 

/a(l)b(l) (-V-  )dV.  =  aK  where 
-Ld  X  s 

aK  =  -2af2e~w(l+2w(e‘6R-(l-6R))/62R2)  (C.6) 

S  * 

Let  B2  -  (1  -  0R  -  e"eR)/62R2 


Then  aK  =  -2af2e"W(l  -  2wB„) 
s  l 

fa2(l)b2(2)V  dV  dV,  =  aJ  '  where 
'  12  1  2  s 

aJg'  =  -aw3Y(e-2w(2w2+2w+l)/4w3  -  e~2w/6w  +  e_2wY(2w+l) /2v3 
-Y((a+.56)R+l)e“(a+‘5B)R/wR2(a+.5e)2 


x(e(a-.56)R_e-(a-.50)R)/(a_<5B)R 


•2fe_(u+-5B)R/wR(a+.56)t(e(a",5B)R-e~(Ct"-5B)R) 


/2(a-.5B)2R2  -(e(a_'56)R+e"(a"*5B)R)/2(a-.5B)R]) 

=  af2e~BR(-w2(4f-l)/3-w(2f2-l)/2  +  (2f 2+2f+l) /4) 

—  RR  ? 

+afve  ptv(.5+w+2w  /3)Y1+2fw(l+2w)Y2 

+8f2w((h+2w)Y3+2wY4))  (C . 7) 

/a(l)b(l)a(2)b(2)V12dV1dV2  =  aKs'  (C.8) 

This  integral  docs  not  have  a  closed  analytic  expression.  The 
procedure  for  numerically  evaluating  this  integral  is  given 


in  appendix  C 
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Procedure  for  evaluating  K' 


aKg’  =  /ls(a1,rla)ls(a2,rlb) (exp(-6r12)/r12)ls(a3>r2a)ls(a4,r2b) 


x  dv^dv2 


ls(a,r)  = 


(a^/ir)1^2  exp(-ar) 


Expand  all  functions  in  spherical  harmonics  about  center  a  using 


(otrb)n  exp(-arb)  =  I  (2*.+lKin(Ctra,aR)Pt(cosea) 

0 


3/2  2 

aKg'  =  BUa^a^)  /tt  )/  expC-o^r^)  exp(-ct,,rlb)  x 


E  <2'-+l)r.,o(Brla,Br2a)Pj(cos0la)  expC-c^r^)  x 


l  ^  ‘2*2+1>^2l(04t2a'“4R)Pt2<c<,se2a)  d“2r2  dr2d''l 

oo 

AB((a1a2a3a4)3/2/u2)  l  (2JI+1)  /expt-c^r^)  exp(-a2rlb)  x 


ttO<6rl.>Sr2«)  ‘‘xP(-V2aKn<V2a’a4R)  x 


P2(cos6la)r2  dr2dvl 


X 


4B((a1a2a3a4)J/^/iT)  E  (2fc+l)  /  expHx^r^) 

2^=0 


£  f0(2V1K£1l(°2rla*a2R)  PJl1(cos9la^Jl0(erla,er2a)  x 
eXp("a3r2aKU(a4r2a'a4R)P£(coseia)r22dr2dWlrl2drl  " 


16B(al02a3a4)  ^(21+1)/  x 


C£0(9rla,9r2a)eXp(“a3r2a^U^a4r2a,a4R^r2  dr2rl  drl 


f’£0^tl*t2^  =  i£^t<^k£^t>^ 

Hl^l’^  *  t>i£^t<^k£-l^t<^  “  t<i£+l^t<^k£^t>^ 

where  i^Cx)  and  k^vx)  are  spherical  modified  Bessel  functions  for 
which  the  simplest  are 

1q(x)  =  sinh  x/x;  i  ^(x)  ■»  cosh  x/x 
kQ(x)  =  k^U)  =  e  x/x 
and  the  recursion  relations  are 

(x)  *=  i£+1(x)  +  ( (22.+1)  /x)  i^(x) 

k£+l^x^  “  k£-l^x^  +  ((2^+l)/x)k^(x) 

The  subroutines  SSIBFN  and  SSKBFN  compute  the  scaled  spherical 


Bessel  Functions 


ig5 (x)  =  ((2£+l)!!/x  )  exp(-x)±£(x) 


k£S(x)  =  (x*'+1/(2£-l)  !  !  )exp(x)k£<x) 


where  (2£+l) ! !  =  1*3* . . .  •  (22.+1)  is  the  "odd  number  factorial", 
In  terms  of  these  scaled  functions 


(2J+DK,  -  (r<l/r>U1)lIs(6r<)kt8(Br>)  * 


exp(-6(r>-r<)) 


CU(a4r2a’a4R)  =  a4  exp(-a4<r>-r<))  x 


[((l/(2M)(2W))(r<V>Ul)i/(a4r<)kws(a4r  ) 


-((l/(2£+l)(2£+2))(r<U2/r>S'+1)i£+1s(a4r<)k£S(a4r>)] 


In  order  to  avoid  negative  indices,  use  is  made  of  kg  =  k  ^  so 
that  for  the  £.  =  0  case 


C01(V2a*a4R)  =  (i0S(a4r<)a4r>_ilS(a4r<)a4r<  /3r>}  X 


exp(_  4(r>~r<)) 
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100  T=3C0\TI=.01*T-1\  REM  . CALCULATE  EL2  DATA 

10*  FOR  C1I=0  TO  4\  FOR  CI=4  TO  8\  GOSUB  102\  NEXT  CI\  NEXT  C1I\  RETURN 
102  DIM  LC1 ( 4 ) , LC ( 8 ) , PF ( 3 ) \  GOSUB  120\  GOSUB  252\  GOSUB  254 
110  GOSUB  200 

112  GOSUB  300\BL=B\  GOSUB  260 

114  IF  ABS < ( BL-B ) /B ) > . 005  THEN  112 

116  GOSUB  130\  GOSUB  400\  GOSUB  140\  RETURN 

120  PRINT  '  T;  LC1;  LC;  LC2;  C1I;  Cl;'? 

121  PRINT  C1I,CI 

122  FOR  1=0  TO  4\LC1 ( I )=15 . 5+. 5*1 \  NEXT  I 

123  FOR  1=0  TO  3\LC( I >=LC1 ( Cl I ) — . 5*1 \  NEXT  I\LC(4)=10 

124  FOR  1=5  TO  7\LC(I >=LC1 (C1I >-4+. 5*I\  NEXT  I 

125  LC<  8 )=LC1 (Cl I)*. 9999999 

126  LC1=LC1(C1I)\LC=LC(CI)\C=10ALC\C1=10ALC1 

1 27  C2=C1-C\LC2=L0G1 0 (1+C2 ) \CT=C1 +C2\R=C1 /CT\PF < 0 )=R*R 

128  PF( 1 ) ,PF< 3 )  =  2*R*( 1-R ) \PF( 2 )  =  ( 1-R >*( 1-R ) \  RETURN 

130  DIM  #2 , XFNPB <  4 , 4 , 4 , 8 )  \  OPEN  'XFNPB'  AS  FILE  #2,  VIRTUAL 

131  XFNPB(0 ,TI ,C1I ,CI )=XF\XFNPB< 1 ,TI ,C1I ,CI >=EF 

1 32  XFNPB ( 2 , TI , Cl I , Cl ) =LN\XFNPB <  3 , TI , Cl I , Cl >  =LP 

133  XFNPB (4,TI,C1I,CI)=B\  CLOSE  \  RETURN 

140  DIM  #3,EVX<3,4,4,8,33) 

141  OPEN  'EV'  AS  FILE  #3,  VIRTUAL  \K=8000\  FOR  Q=0  TO  3 

142  FOR  1=0  TO  32\EVX(Q,TI ,C1I ,CI , I )=K*EL(Qf I ) 

143  NEXT  I\  NEXT  Q\  CLOSE  \  RETURN 

138  OPEN  'LP:'  AS  FILE  #1\  PRINT  #1 ,DATE*<0X) ,TIME*(0) 

199  CLOSE  \  RETURN 


DOC  REM  ...CALC.  PRELIM.  EF,  LN ,  LP  &  B.  (Host  =  GaAs ) 

202  DIM  E<  2 '  \  GOS'JB  232\  IF  FF<>51  THEN  GOSUB  246 

204  KT=. 0661 7*T/ED\XK=2*EXP( 1 /XT )\N0=4. 8293*1 0A15\B0= . 0209983*N0 
206  NC=N0*MC*'ME*T>A1 .5\NV=N0*MV*<MH*T)Ai.5\  IF  LC>9  THEN  212 

203  ZF= . 5  * ( EG+KT*LOG ( NV/NC ) > \XF=EXP < EF/KT ) \LN=L0G1 0 ( NC ) +LQG10 (XF > 

210  3=10A-6\LP=LOG10 ( NV>+. 4343* <  EG-EF ) /KT\  GCSUB  250\  RETURN 

2! 2  CC=C/NC\  IF  CC<.0001  THEN  XF=CC/ < l-CC/4 ) V  GOTO  216 
214  ?M= . 5* ( l-CC/4  t /SQR( CC*XK ) \XF=SQR < CC/XK ) * < S0R ( TM*TM+ 1 ) -TM  > 

216  FOR  M=0  TO  2\  GOSUB  224\E(M>=XF\  NEXT  M\  GOSUB  244 
216  GOSUB  240\LN  =  L0G10(NO+L0G10(XF)+L0G10(Sl  )\  GOSUB  242 
220  B=SQR i B0*MC*MEA 1 . 5*SOR <T>*S2/KE) *RB\ EF=KT*LOG (XF ) 

222  LP= LOGIC <  NV  H- . 4343*  < EG-EF > /KT\  GOSUB  250\  RETURN 

224  GOSUB  240\Y=EG/KT\  GOSUB  280 

226  '<?= '  MV/MC  >  *  <  MH/ME )  A 1  .  5*EXP< Y  )  /  <  XF*S1 ) 

223  XF=CC/S1/(1+XK*XF>+KP\  RETURN 

230  ED=44\RB=20*10A-6\EG=- ( 1 160- . 7Q2*TA2/ (T+1108 ) ) /ED 

231  m:=£\MV=1\KE=11. 8\ME= . 327\MH= . 559\  RETURN 

232  ED=30\EA=28\R8=40*10A-8\£G=-'1522-.5B*TA2/(T+300) )/ED 
223  MC=1\MV=1\KE=13.1\ME=.067\MH=.473\  RETURN 

236  ED-3 2\RB=37*10a-8\EG=-(741-.456*Ta2/(T+210> >/EB 

23"  MC=5\M'.'=1  \KE=1 6\ME= .  221\MH=.  31  \  RETURN 

242  S 1 =  0 \  FOR  0=0  TO  FF\S1=S1+NX ( J ) / (XF+EX ( J ) ) \  NEXT  J\  RETURN 

242  S2=D\  FOP  J=Q  TO  FF\D=XF+EX< J)\S2=S2+NX( J)*EX< J)/(D*D> 

24  3  NEXT  JV 32=XF*S2\  RETURN 

244  TM=2*Et 1 >  — E ( 0 ) — E ( 2 > \  IF  TM=0  THEN  XF*E(2)\  RETURN 

245  XF=E<  2 ) +ABS< E< 1 )-E<  2  3 ) * 2/TM\  RETURN 

246  FF,I  =  51\D=. 0 2 \  DIM  EX< I ) , NX< I ) \A=2*D/3/SQR(PI ) \  FOR  K=1  TO  26 
24"  I=2*K-2\X=K*K*D\EX(I >=EXP<X)VNX<I>=SQR(X>*4*X*A\X=KA(K+1)*D 

249  EX'  1  +  1 ' =EXP'X ) \NX( 1+1 )=SQR(X>*( 2*K+1 >*A\  NEXT  K\  RETURN 

250  OPEN  'LP:'  AS  FILE  #1\  PRINT  # 1 , LC1 ; LC ; B , EF , LN , LP\  CLOSE 

251  PRINT  LC1 ; LC , B , EF , LN , LP\  RETURN 

252  OPEN  'LP:'  AS  FILE  #1\  PRINT  #1 f 'Temp' ; 'LogCl ' ; 'LogC' ; 'LogC2' 

253  PRINT  f 1 ,T; LC1 ; LC; LC2\  CLOSE  #1\  PRINT  T; LC1 ; LC; LC2\  RETURN 

254  OPEN  'LP:'  AS  FILE  #1\  PRINT  #1,'LC1  'LC ','B' , 'EF  \'LN  VLP 

255  CLOSE  #1\  PRINT  'LC1  ' ; 'LC' , 'B' , 'EF' , 'LN' , 'LP'\  RETURN 

260  REM  _ RECALCULATE  EF,  LN,  LP  &  B,  USING  EL(Q,I) 

262  CB=C*REA3\CR=CB*4*PI/3\  IF  SF<>51  THEN  GOSUB  246 
270  E(0)=XF\  FOR  M=1  TO  2\S=0\AC=1\  GOSUB  274\  GOSUB  224 
272  GOSUB  290\E(M)=XF\  NEXT  M\  GOSUB  244\  GOSUB  218V  RETURN 
274  FOR  1=1  TO  32\R=.25*I\  GOSUB  282\  NEXT  I\I=33 
280  A=AC\EL<  0,33)  ,EL(  1 , 33 )=. 5*(EL( 0 ,32 )+EL(l ,  32 ) ) \  GOTO  284 
262  AL=AC\Y=-CR*RA3\  GOSUB  288\AC=EXP(Y)\A=AL-AC 

294  Y=EL(0,I)/KT\  GOSUB  20B\Y1=.5*EXP(Y)/XF 

295  Y=EL(1,I>/XT\  GOSUB  28B\Y2=.5*EXP(-Y)*XF 
286  S=3+A* ( . 5+Y1 ' / ( 1 +Y1 +Y2 ) \  RETURN 

2S9  IF  Y<-40  THEN  Y=-40\  RETURN  ELSE  RETURN 

290  XF=CC*S/S1+KP\  RETURN 


300  REM  _ ENERGY  EL(Q)  VS  R  (SCRN.  POT.,  VAR.)  1/86 

302  IF  KFOl  THEN  GOSUB  364 

304  FOR  0=0  TO  l\A=<2-Q*5/16)*<l-BA2/8)\AC=A+.06 
306  FOR  1=0  TO  32\R=.25*I\  GOSUB  330 

308  AT=A\  FOR  M=0  TO  2\A=AT+. 08*(M-1 ) \  GOSUB  336\E(M)=E\  NEXT  M 

310  GOSUB  31 8\  IF  ABS( A-AT) > . 02  THEN  308  ELSE  GOSUB  336\AL=AC\AC=A 

312  IF  0=0  THEN  EL(0,I)=E\  IF  EL<0,I)>0  THEN  EL(0,I>=0 

314  IF  Q=1  THEN  EL( 1 ,1 )=E-EL< 0 , 1 ) \  IF  EL(1,I)>0  THEN  EL(1,I>=0 

316  A=2*AC-AL\  GOSUB  322\  NEXT  I\  NEXT  Q\  RETURN 

318  DX=E < 0 ) -2*E ( 1 ) +E { 2 > \  IF  DX<.01  THEN  DX=.01 

320  A=AT+. D4*(E(0)-E<2) )/DX 

322  IF  A <.01  THEN  A=.01\  RETURN  ELSE  RETURN 

324  E=,Q+1)*(AA-4*A*HH)+Q*.5*A*HH*(.5+H+HH>\  RETURN 

326  E= <  AA* ( 1 +SK ) +A* ( 2*K-2*HH+ J ) ) / < 1+S ) \  RETURN 

328  E=<2*AA*<1+S*SK)+A*(2*J-4*HH+JP+4*S*K+KP) )/<l+S*S>\  RETURN 

330  U=B*R\BE=EXP ( -U )  \  IF  U>.2  THEN  B2=<1-BE-U)/(U*U)\  RETURN 

332  E2=-(360-U*(120-U*<30-U*(6-U) ) ) )/720\  RETURN 

336  AA=A*A\H=A/(A+.5*B)\HH=H*H\W=A*R\WW=W*W\  IF  R=0  GOTO  324 

338  W3=WW/3\AX=EXP<-W)\Y=2*W-U\YY=Y*Y 

34C  S=AX*a+W+W3)\SK=S-2*AX*W3\  IF  YY>.25  THEN  344 

342  Y4=-<1680-336*Y+56*YY-8*Y*YY+YY*YY>/40320\  GOTO  346 

344  Y4=( 1-EXP(-Y)-Y+. 5*YY-Y*YY/6 )/YY/YY 

346  Y3=Y*Y4+1 /6\Y2=Y*Y3- . 5\Y1=Y*Y2+1 \K=-2*HH*AX* <1-2*W*B2 ) 

348  J=-2*H*BE*(2*W*<2*H*Y2+Y1)+H)\  IF  Q=0  GOTO  326 

350  TM=< .5+W+2AW3)*Yl+2*HA,W+6*W3)*Y2+8*HH*W*( (H+2*W)*Y3+2*W*Y4 ) 

352  0P=BE* ( W*H*TM+HHA ( . 5* ( HH+H+ . 5 ) -W* { HH- . 5 ) -W3* ( 4*H-1 ) > ) 

354  3=ABS(5*B/A)\X=INT(G)\D=G-X\  GOSUB  360\SM=TM*(1-D)\X=X+1 

356  GOSUB  360\SM=SM+TM*D\KP=.5*HH*< .5+H+HH)*AX*AX*(l+W*SM)\  GOTO  32B 

360  P=X\  IF  P>24  THEN  P=24 

362  TM=C1(P)+W*(C2(P)+W*<C3(P)-WAC4(P)))\  RETURN 

364  DIM  EL <3,33) ,LE(2) ,A(2) 

365  DIM  Cl < 24 ) ,C2(24) ,C3<24) , C4 < 24 ) \  RESTORE  367 

365  FOR  P=0  TO  24\  READ  Cl <P) ,C2(P) ,C3(P) ,C4 (P)\  NEXT  P\KF=1\  RETURN 

367  DATA  2 . 00243, 1 . 58709 55446,- . 01856 , 2 . 00196, 1 . 56891 51171 00428 

368  DATA  2 . 00011 , 1 . 55993 46202 ,. 01720 , 1 . 99791 , 1 . 55498 ,. 41181 ,. 02403 

369  DATA  1 . 99577 , 1 . 55152, . 36420 ,. 02712, 1 . 99391 , 1 . 54822 32063 02795 

370  DATA  1.99240,1.54452, .28158, .02743, 1.991 24, 1.54018,. 24706,. 02616 

371  DATA  1.99038,1.53523, .21676, .02449,1.98978,1.52975, .19029, .02265 

372  DATA  1 . 98939, 1 . 52385 ,. 16723, . 02078, 1 . 98915 , 1 . 51768, . 14716, . 01895 

373  DATA  1 . 98904 , 1 . 51133 ,. 12970 ,. 01721 , 1 . 98900, 1 . 50492, . 11449 ,. 01559 

374  DATA  1 . 98902, 1 . 49854 ,. 10123 ,. 01410, 1 . 98906 ,1 . 49224 ,. 08964 ,. 01273 

375  DATA  1 . 98912 , 1 . 48609 ,. 07949 ,. 01149 , 1 . 98917 , 1 . 48013 ,. 07057 ,. 01036 

376  DATA  1.98921,1.47438, .06273, .00933,1.98923,1.46887, .05580,. 00840 

377  DATA  1 . 98924 , 1 . 46361 ,. 04966 ,. 00755 , 1 . 98922 , 1 . 45859, . 04422 ,. 00678 

378  DATA  1.90918,1.45383, .03936, .00609,1.98912,1.44931,  .03502, .00546 

379  DATA  1 . 98905 , 1 . 44504 ,. 031 1 4 ,. 00489 


400  REM  _ ENERGY  El  &  E2  VS  R  (COMPENSATED)  8/84 

402  IF  KFOl  THEN  GOSUB  364 
404  DA=. 04\DL=. 04 

406  FOR  Q=2  TO  3\A=1-BA2/8\L=0\E=-1 

408  FOR  1=32  TO  0  STEP  -1\  IF  E>=0  THEN  E,A,L=0\  GOTO  414 
410  R=.25*I\  GOSUB  420 

414  EL ( 0 , I ) =E- ( Q-2 ) *EL ( Q-l , I ) \  IF  EL(Q,I)>0  THEN  EL(Q,I)=0 
416  NEXT  I\  NEXT  Q\  RETURN 

420  LT=L\  FOR  N=0  TO  2\L=LT+(N-1 )*DL\  GOSUB  430\LE(N)=E\A(N)=A\  NEXT  N 
422  D=LE ( 0 ) -2*LE ( 1 ) +LE < 2 > \  IF  D<.001  THEN  D=.001 
424  L=LT+.5*DL*(LE(0)-LE(2) )/D 
426  IF  ABS(LT-L) >DL  THEN  420 

428  A=A(1)+.5*(L-LT)*<A(2)-A(0) )/DL\  GOSUB  440\  RETURN 

430  AT=A\  FOR  M=0  TO  2\A=AT-(M-1 )*DA\  GOSUB  440\E(M)=E\  NEXT  M 

432  D=E(0)-2*E<1 >+E<2>\  IF  D<.001  THEN  D=.001 

434  A=AT+.5*DA*(E(2)-E(0> )/D\  IF  A<DA+.001  THEN  A=DA+.001 

436  IF  ABS(AT-A) >DA  THEN  430  ELSE  GOSUB  440\  RETURN 

440  AA=A*A\H=A/(A+.5*B)\HH=H*H\X=ABS(L>\  GOSUB  460\  GOSUB  466\J1=J 

442  X=ABS(L+R)\  GOSUB  460\  GOSUB  466\J2=J\  IF  Q=2  THEN  456 

444  X=ABS(2*L+R)\  GOSUB  460\W3=W*W/3\  IF  YY>.25  THEN  44B 

446  Y4=-<1680-Y*(336-Y*<56-Y*(8-Y) ) ) ) /40320\Y3=Y*Y4+l/6\  GOTO  450 

448  Y4= ( 1-EXP ( -Y ) -Y+ . 5*YY-Y*YY/6 ) /YY/YY\Y3=Y*Y4+1 /6 

450  TM=( .5+W+2*W3)*Yl+2*H*(W+6*W3)*Y2+8*HH*W*< <H+2*W)*Y3+Z*W*Y4 ) 

452  JP=BE* < W*H*TM+HH* ( . 5* ( HH+H+ . 5 ) -W* ( HH- . 5 ) -W3* ( 4*H-1 ) >  > 

454  E=2*AA+A*(2*J1-2*J2-JP>\  RETURN 

456  E=AA+A*( Jl-J2)\  RETURN 

460  U=B*X\BE=EXP ( -U ) \W=AAX\Y=2*W-U\YY=Y*Y 

462  IF  YY> . 25  THEN  Y2= (1-EXP (-Y)-Y)/YY\Y1=Y*Y2+1\  RETURN 

464  Y2=-(36Q-Y*(120-Y*(30-Y*(6-Y) ) ) )/720\Yl=Y*Y2+l\  RETURN 

466  J=-2*H*BE*(2*W*(2*H*Y2+Y1)+H)\  RETURN 

490  PRINT  L,A,E(0) ,E(1 ) ,E(2) ,N\  RETURN 

491  PRINT  L,A, ,E, ,R\  RETURN 

432  PRINT  N , , LE ( 0 ) , LE ( 1 ) , LE ( 2 ) \  RETURN 
493  PRINT  L,A, ,E, ,R\  RETURN 


fc  *  • 


500  IF  QX<>3  THEN  GOSUB  560\  REM  _ CALC.  DS2  DATA 

502  T=300\TI=.01*T-1\  FOR  C1I=0  TO  4\  FOR  CI=4  TO  8\  GOSUB  504 

503  NEXT  CI\  NEXT  Cl I \  RETURN 

5C4  FOR  1=0  TO  4\LC1 < I >=15 . 5+. 5*I\  NEXT  I 

505  FOR  1=0  TO  3\LC(I >=LC1 (C1I >-. 5*I\  NEXT  I\LC<4>=10 

506  FOR  1=5  TO  7\LC(I)=LCl<ClI)-4+.5*I\  NEXT  I 

507  LC(8)=LC1(C1I>*. 9999999 

508  LC1=LC1 (C1I >\LC=LC<CI ) \C=10ALC\C1=1CALC1 

509  C2=C1-C\LC2=LOG10(1+C2 ) \CT=C1+C2\R=C1/CT\PF< 0 >=R*R 

510  PF ( 1 ) , PF ( 3 ) =2*R* ( 1-R ) \PF ( 2  > = ( 1 -R ) * ( 1-R ) 

520  GOSUB  530\  GOSUB  800\  GOSUB  700\  GOSUB  810\  GOSUB  710 

521  GOSUB  660\  PRINTSCREEN  \  GOSUB  600\  GOSUB  B20\  GOSUB  720 

522  GOSUB  830\  GOSUB  730\  GOSUB  605\  GOSUB  735 

523  PRINTSCREEN  \  GOSUB  61 0\  GOSUB  840\  GOSUB  740 

524  GOSUB  790\  PRINTSCREEN  \  GOSUB  597\  RETURN 

530  OPEN  'XFNPB'  FOR  INPUT  AS  FILE  *2,  VIRTUAL 

531  OPEN  'EV'  FOR  INPUT  AS  FILE  #3,  VIRTUAL 

532  EF=XFNPB < 1 , TI , Cl I , Cl ) \ LN=XFNPB <  2 ,  TI ,  Cl  I ,  Cl ) 

533  LP=XFNPB ( 3 , TI , Cl I , Cl ) \B=XFNPB  <  4 , TI , Cl I , Cl ) \K= . 0001 25 

534  XF=XFNPB(0,TI,C1I,CI)\  FOR  Q=0  TO  3\  FOR  1=0  TO  32 

535  EL ( Q , I ) =K*EVX ( Q , TI , Cl I , Cl , I ) \  NEXT  I\  NEXT  Q\  CLOSE 

536  ES=EL( 1 , 32)\  IF  ES<-1  THEN  ES=-1 

537  FOR  Q=0  TO  3\EL(Q, 33)=ES\  NEXT  Q\EF<2> ,EF<3)=. 8 

538  EF<  0 ) ,EF(1 )=1\RF(0 ) ,RF(1 )=1\RF(2) ,RF(3)=1 .25\  GOTO  542 

540  ED=44\RB=20*10A-8\EG=-(1160-.702*TA2/mil08>)/ED 

541  MC=G\MV=1 \KE=11 . 8\ME= . 327\MH= . 559 

542  ED=30\EA=28\RB=40*10A-8\EG=- ( 1522- . 58*TAZ/ (T+300 ) ) /ED 

543  MC=1\MV=1 \ KE=1 3 . 1 \ME= . 067\MH= .473 

544  KT=.08617*T/ED 

545  CB=CT*RBA3\CR=CB*4*PI / 3\M=2/ 3\PM= ( CR/M ) AM*EXP (M ) 

546  FOR  1=0  TO  32\R=. 25*I\Y=-CR*RA3\  GOSUB  648 

547  P=PM*RA2*EXP<Y>\P{I)=P\R<I)=R\RH(I)=RF(Z>*R\  NEXT  I 
543  RM= ( 2/3/CR > A ( 1/3 ) \FGX=100*TI+10*C1I+CI+100 

550  N0=4 . 8293*1 0A15\NC=N0*MC*( ME*T ) A1 . 5\NV=N0*MV*<MH*T) A1 . 5 

551  LNCO=LOG1 0  <  NC  > -1 . 5*L0G1 0 ( KT ) \ LNVO=LOG1 0  <  NV ) -1 . 5*L0G1 0 ( KT ) 

552  RETURN 

560  DIM  EL(3,33) ,P(32) ,PF(3) ,EF(3) ,RF(3) ,R(3Z) ,E(200) ,X<3) 

561  DIM  L<( 3)\  DIM  #2 , XFNPB ( 4 , 4 , 4 , 8 ) \  DIM  #3, EVX(3, 4,4,8,33) 

562  FOR  1=0  TO  200\E(I ) = . 02*I\  NEXT  I\QX=3\IX=200\RX=32\LX=13 

563  DIM  LC1 < 4  > , LC <  8 ) , Y ( RX > , S ( IX ) , RH ( RX ) , L ( LX ) 

564  DIM  #4 , SFX (1,4,4,8,13)\  DIM  f 5 ,SDX< 3 , 4 , 4 , 8 , 13 > 

565  DIM  #6 , SAX (3,4,4,8,13)\  RETURN 

570  OPEN  'SF'  AS  FILE  #4,  VIRTUAL  \K=1000\  FOR  1=0  TO  LX 

571  SF*<Q,TIrClI,CI,I)=K*L(I>\  NEXT  I\  CLOSE  \  RETURN 

572  OPEN  'SD'  AS  FILE  #5,  VIRTUAL  \K=1000\  FOR  1=0  TO  LX 

573  SDX ( 0 , TI , Cl I , Cl , I ) =K*L < I > \  NEXT  IN  CLOSE  \  RETURN 

574  OPEN  'SA'  AS  FILE  #6,  VIRTUAL  \K=1000\  FOR  1=0  TO  LX 

575  SAX (0,TI,C1I,CI,I) =K*L ( I ) \  NEXT  I\  CLOSE  \  RETURN 

597  OPEN  'LP:'  AS  FILE  #1\  PRINT  #1,CHR$(12)\  CLOSE  \  RETURN 

598  OPEN  'LP:'  AS  FILE  #1\  PRINT  #1 , DATE$ ( OX ) , TIME* ( OX ) 

599  CLOSE  \  RETURN 

281 


V'V  v-  .  ^  V 


600  DIM  DS(QX,IX>\  REM  ...  DEN.  OF  <DD>  ELECTRON  STATES 

601  SW=3\F=0\  FOR  0=0  TO  1\  GOSUB  620\  NEXT  Q\  REM  <— ) 

602  F=l\  FOR  Q=2  TO  3\  GOSUB  620\  NEXT  Q\  GOTO  616\  REM  <-+> 

605  DIM  DS<QX, IX) \  REM  ...  DEN.  OF  FILLED  (DD)  STATES 

606  SW=1\F=0\  FOR  Q=0  TO  1\  GOSUB  620\  NEXT  Q\  GOTO  616 
610  DIM  DS(QX,IX)\  REM  ...  DEN.  OF  HOLE  STATES 

612  SW=3\F=2\  FOR  Q=0  TO  1\  GOSUB  620\  NEXT  0\  REM  {++) 

614  F=3\  FOR  Q=2  TO  3\  GOSUB  620\  NEXT  0\  REM  (+-) 

616  GOSUB  652\  RETURN 

620  A2=PF ( F ) *CT\AC=0\E2=-EF ( F ) *EL ( Q , 0 ) \ EC=I NT ( 50*E2 ) 

622  FOR  1=0  TO  33\R=RF(F)*. 25*(I+.5)\TM=CR*R*R*R\E1=E2 

624  E2=-EF(F)*EL(Q, I ) \DX=ABS(E2-E1 ) \  IF  DX<.02  THEN  DX=.02 

626  A1=A2\  IF  TM<25  AND  I<33  THEN  A2=CT*PF(F)*EXP(-TM>  ELSE  A2=0 

626  AL=AC\AC= ( A1-A2 ) /DX\  IF  SW=3  THEN  636 

630  Y=EF(F)*EL(0,I)/KT\  GOSUB  648\Y1=. 5*EXP(Y)/XF 

632  Y=EF(F)*EL(1,I)/KT\  GOSUB  648\Y2=. 5*EXP<-Y)*XF 

634  IF  Q=0  THEN  AC=AC*U+Y2)/(1+Y1+Y2)  ELSE  AC=AC*Y2/ < 1+Y1+Y2 > 

636  EL=EC\EC=INT( 50*E2 ) \D=EC-EL\  IF  1=0  THEN  646 

638  IF  D=0  THEN  E=EL\A=AC\  GOSUB  650\  GOTO  646 

640  M= ( AC-AL ) / ABS < D ) \  IF  D<0  THEN  D=-D\  GOTO  644 

642  FOR  X=1  TO  D\E=EL+X\A=AL+X*M\  GOSUB  650\  NEXT  X\  GOTO  646 

644  FOR  X=1  TO  D\E=EL-X\A=AL+X*M\  GOSUB  650\  NEXT  X 

646  NEXT  I\  RETURN 

648  IF  Y<-40  THEN  Y=-40\  RETURN  ELSE  RETURN 

65C  DS  <  Q , E ) =A+DS < Q , E ) \  RETURN 

652  L=.25\M=.5\N=.25\  FOR  Q=0  TO  SW 

654  Y=DS(Q,0)\Z=DS(Q,1 )\  FOR  1=2  TO  IX\X=Y\ Y=Z\Z=DS ( Q , I ) 

656  DS(Q,I-1)=L*X+M*Y+N*Z\  NEXT  I\  NEXT  Q\  RETURN 

660  SET  LINE  STYLE  4\D=.125*B\  FOR  1=10  TO  32 

661  PLOT  <.25*I,-32*EXP<-D*I)/I>,\  NEXT  I\  SET  LINE  STYLE  1 

662  X=.25*10\Y=-32*EXP(-D*10)/10-.4\L*(3)='Vb'\  GOSUB  728\  RETURN 

670  SET  LINE  STYLE  2\  FOR  1=2  TO  74  STEP  2\  GOSUB  674 

671  PLOT  <~E(I),S),\  NEXT  I\  PLOT  \  SET  LINE  STYLE  1 

672  FOR  1=1  TO  75\  GOSUB  6?3\  PLOT  (-E(I),S)f\  NEXT  I\  PLOT  \  RETURN 

673  E=.02*I\S=A+.5*LOG10(E)\  RETURN 

674  E=. 02*1 \S=A+. 5*L0G1 0(E) -LOG10(1+EXP< (E-EF) /KT> ) 

675  IF  S<14  THEN  S=14\  RETURN  ELSE  RETURN 

700  DIM  Y(RX)\X=RM-.3\  IF  X>7.25  THEN  X=7.25\  REM  . PROBABILITY 

701  L$<  0 )  =  ' AA ' \L$< 1 )  =  '  DA'\L*<2)='  DD' 

702  FOR  Q=0  TO  2\  FOR  1=0  TO  RX\Y( I )=Y(I )+PF(2-Q)*P(I >\  NEXT  I 

703  Y=Y ( 4*X+2 ) — . 24 \  GOSUB  717\  GOSUB  728\  NEXT  Q\X=4\Y=1 

704  Li<2)='(a)  Impurity  Spacing  Distribution' \  GOSUB  728\  RETURN 


710  L$(0)='E1'\L*(1)='E2'\L$(2)='E3\  REM  . DONOR  LEVELS 

711  L$( 3)='E4 '\X(0)=3\X( 1 )=1\X(2)=8\X(3 )=5 

712  FOR  0=0  TO  3\  FOR  1=0  TO  RX\Y(I )=EL(Q,I >\  NEXT  I 

713  J=X(G>\X=R(J)\Y=Y(J>-.4\  GOSUB  717\  GOSUB  728\  NEXT  Q 

714  X=3\Y=-3. 8\L$(3)=' (b)  Donor  Levels  <Ryd>'\  GOSUB  728\  RETURN 

715  SET  LINE  STYLE  2\  FOR  1=0  TO  200  STEP  2\  PLOT  <E(I>,S<I>>, 

716  NEXT  I\  PLOT  \  SET  LINE  STYLE  1\  RETURN 

717  FOR  1=0  TO  32\  PLOT  <R< I ) , Y( I ) > , \  NEXT  I\  PLOT  \  RETURN 

718  FOR  1=0  TO  32\  PLOT  (RH(I ) ,Y(I ) ) ,\  NEXT  I\  PLOT  \  RETURN 

720  L$(0)='HI ' \L$(1 )='H2' \L$( 2 )='H3'\  REM  . ACCEPTOR  LEVELS 

721  L$t  3 )  =  'H4 '  \  FOR  Q=0  TO  3\  FOR  1=0  TO  RX\Y(I )=-EF(2)*EL(Qf I ) 

722  NEXT  I \ J=X ( 0 ) \X=RH <  J ) \ Y=Y ( J ) \  GOSUB  718\  GOSUB  728\  NEXT  Q 

723  X=3\Y=2.8\L*(3)='(c>  Acceptor  Levels  (Ryd)'\  GOSUB  728\  RETURN 

725  IF  <Q=0)  OR  (0=2)  THEN  PLOT  CURVE  (E( ) ,S( ) ,201,0)\  RETURN 

726  SET  LINE  STYLE  4\  FOR  1=0  TO  200  STEP  2\  PLOT  <E(I),S(I)), 

727  NEXT  I\  PLOT  \  SET  LINE  STYLE  1\  RETURN 

728  SET  POSITION  <X,Y)\  GRAPHIC  PRINT  L$(Q)\  RETURN 

730  L*(0)='E1'\L*(1 >='  E2'\  REM  . DONOR  STATES 

731  L*(2)='E3'\L$<3)='E4'\  FOR  0=0  TO  SW\  GOSUB  750 

732  GOSUB  725\  GOSUB  728\  GOSUB  572\  NEXT  Q\A=LNC0\  GOSUB  670 

733  X=-1.2\Y=18\L*(SW)='Cond.  Band'\  GOSUB  728\X=2\Y=19 . 5 
-734  L*(SW)='(d)  Donor  States  vs  E  (Ryd>'\  GOSUB  728\  RETURN 


735  FOR  Q=0  TO  SW\  GOSUB  750\  REM  . . FILLED  STATES 

736  GOSUB  715\  GOSUB  570\  NEXT  Q\  RETURN 

740  L*(0)='H1'\L*<1)='H2'\L$<2)='H3'\  REM  . ACCEPTOR  STATES 


741  L*  <  3 )  = '  H4'\  GOSUB  745\  FOR  0=0  TO  SW\  GOSUB  750\  GOSUB  725 

742  GOSUB  728\  GOSUB  574\  NEXT  Q\A=LNV0\  GOSUB  672\  GOSUB  745 

743  X=. 3\Y=18\L$(SW)= 'Valence  Band'\  GOSUB  728\X=-4.5\Y=19.5 

744  L*  <  SK  >  = '  ( e )  Acceptor  States  vs  E-EV  <Ryd)'\  GOSUB  72B\  RETURN 

745  FOR  1=0  TO  IX\E< I )=-E< I ) \  NEXT  IN  RETURN 

750  Y=14\  FOR  1=0  TO  IX\S( I ) =LOG10 ( 100+DS(Q , I ) ) 

751  IF  S < I ) <14  THEN  S(I)=14  ELSE  IF  Y<S(I)  THEN  Y=S(I)\IM=I 

752  NEXT  I\X=E(IM)\  DIM  L(LX)\L(0)=Y\L(7)=X\L2=INT(Y) 

753  L=6\  FOR  I=IM  TO  0  STEP  -1 

754  IF  S< I ) >L2-6+L  THEN  755  ELSE  L(L)=E(I >\L=L-1\  IF  L>=1  THEN  754 

755  NEXT  I\L=8\  FOR  I=IM+1  TO  IX 

756  IF  S(I ) >L2+8-L  THEN  757  ELSE  L(L)=E(I )\L=L+1\  IF  L<=13  THEN  756 

757  NEXT  I\  RETURN 

790  SET  CHARACTER  SIZE  .12,. 6\  SET  CHARACTER  SPACING  .12,0 

791  SET  POSITION  (-6 . 65 , 12 . 5 ) \  GRAPHIC  PRINT  USING  'Fig.  A  #*#.', FGX; 

792  FOR  Y=1  TO  5\  GOSUB  793\  SET  POSITION  (-6 . 25 , 12 . 5-. 6*Y) \  NEXT  Y\  RETURN 

793  ON  Y  GOSUB  794, 795, 796, 797, 798\  GRAPHIC  PRINT  L*\  RETURN 

794  L*='  (a)  Probability,  AA  +  DA  +  DD  ,  of  Pairs  vs  R'\  RETURN 

795  L$='(b)  Donor  Energy  E1;(DD)-,  E2;(DD-)-,  E3; ( DA > - ,  E4(DA+>-  vs  R'\  RETURN 

796  L$='(c>  Acceptor  En.  H1;<AA>  +  ,  HZ; < AA+  >  +  r  H3;(AD)  +  ,  H4(AD-)+  vs  R'\  RETURN 

797  L$='<d)  Densities  of  Donor  States,  El,  E2,  E3,  &  E4,  vs  E'\  RETURN 
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800  CLEAR  \  SET  VIEWPORT  . 1  , 1 , 0 , . 625\A1=. 03\  REM  . PROBABILITY  AXES 

801  SET  WINDOW  -.2,8,-3f2\  SET  CHARACTER  SIZE  . 1 , . 25\Y1=0\Y2=1 

802  SET  CHARACTER  SPACING  .1,0\  GOSUB  850 

803  SET  POSITION  <-.l,.5>\  GRAPHIC  PRINT  'P'\  RETURN 

810  SET  WINDOW  - . 2 , 8 ,-4 . 2 , 3. 8\Y1=-4\Y2=0\A1=-. 05\  REM  DONOR  ENERGY  AXES 

811  GOSUB  852\  GOSUB  860\  RETURN 

820  CLEAR  \  SET  WINDOW  -. 2 , 8 ,-4 . 6 ,3 . 4\A1= . 05\  REM  ACCEPTOR  ENERGY  AXES 
921  Y1=0\Y2=3 . 6\  GOSUB  850\  GOSUB  B80\  RETURN 

830  SET  WINDOW  -1 . 7,6. 5,13. 1 ,27. 1\  REM . DONOR  DENSITY  OF  STATE  AXES 

831  X1=-1.5\X2=6.5\  GOSUB  870\  GOSUB  845 

832  SET  POSITION  <5.2,14.2>\  GRAPHIC  PRINT  'E  (Ryd)' 

933  SET  POSITION  (0,19.5)\  GRAPHIC  PRINT  'EC'\  RETURN 

840  CLEAR  \  SET  WINDOW  -6. 7 ,1 . 5 , 6 . 1 ,20. 1\  REM  ACCEPTOR  DEN.  OF  ST.  AXES 

841  Xl=-6 . 5\X2=1 . 5\  GOSUB  870 

342  SET  POSITION  <0,19.5>\  GRAPHIC  PRINT  USING  'EV  =  ###.###', EG 
843  SET  POSITION  <-5.95,14.2>\  GRAPHIC  PRINT  'E-EV  (Ryd)'\  RETURN 

845  SET  LINE  STYLE  4\  PLOT  (-EF, 13.4 ) , <-EF,19)\  SET  LINE  STYLE  1 

846  GRAPHIC  PRINT  'EF'\  RETURN 

350  FOR  X=-0  TO  8\  SET  POSITION  <X-.15,-7*A1)\  GRAPHIC  PRINT  X\  NEXT  X 
951  SET  POSITION  (6. 15 ,-9*Al ) \  GRAPHIC  PRINT  'R  (Bohr)'\  PLOT 

852  PLOT  (0,Q),(8,0)\  FOR  X=0  TO  8\  PLOT  (X , 0 ) , (X,3*A1 ) 

853  FOR  Y= . 1  TO  1  STEP  .1\  PLOT  <X+Y,0> , <X+Y,A1 )\  NEXT  Y\  NEXT  X 

854  FOR  X=. 5  TO  B\  PLOT  (X,0) , <X,2*A1 >\  NEXT  X\  PLOT  <0,Y1 ) , (0,Y2) 

855  FOR  Y=Y1  TO  Y2\  PLOT  <0,Y> , < . 1 ,Y)\  NEXT  Y\  FOR  Y=Y1  TO  Y2  STEP  .2 

856  PLOT  <  0 , Y ) , ( . 03 , Y ) \  NEXT  Y 

857  SET  POSITION  <-. 1 , Yl+2. 5 >\  GRAPHIC  PRINT  'E'\  RETURN 

860  SET  POSITION  (6,-2.05)\  GRAPHIC  PRINT  USING  'TEMP.  *  ###  K',T 

861  SET  POSITION  (6,-2. 4)\  GRAPHIC  PRINT  USING  'LOG  Cb  =  ###.#', LC1 

862  SET  POSITION  (6,-2.75)\  GRAPHIC  PRINT  USING  'LOG  C  =###.*', LC 

863  SET  POSITION  (6,-3.1)\  GRAPHIC  PRINT  USING  'LOG  Cx  =  ##*.##*', LC2 

864  SET  POSITION  (6,-3. 45}\  GRAPHIC  PRINT  USING  'Log  n  =  ###.##♦', LN 

365  SET  POSITION  (6,-3. 3>\  GRAPHIC  PRINT  USING  'Log  p  =  ###.#*# ',LP 

866  SET  POSITION  (-.2,-4)\  GRAPHIC  PRINT  -4 

867  SET  POSITION  <-.2,-2>\  GRAPHIC  PRINT  ~2\  RETURN 

B^O  PLOT  (XI, 14) , (X2,14)\  FOR  X=Xl+.5  TO  X2\  PLOT  (X , 14 ) , (X , 14 . 3 ) 

971  NEXT  X\  FOR  X=X1  TO  X2\  PLOT  (X, 14 ) , (X, 14 . Z) \  NEXT  X 
672  FOR  X=X1  TO  X2  STEP  .1\  PLOT  (X,14 ) , (X,14 . 1 )\  NEXT  X 

873  FOR  X=Xl+.5  TO  X2\  SET  POSITION  (X-.l,13.4)\  GRAPHIC  PRINT  -X\  NEXT  X 

874  PLOT  (0,14) , (0,19. 5)\  PLOT  (XI ,14 ) , (XI , 20 > \  FOR  Y=14  TO  19 

875  PLOT  (XI ,Y) , (X1+. 08 , Y ) \  SET  POSITION  <X1-.3,Y>\  GRAPHIC  PRINT  Y 

876  FOR  X=1  TO  4\S=Y+L0G10(2*X)\  PLOT  (XI ,S) , (X1+.04,S)\  NEXT  X\  NEXT  Y 

877  SET  POSITION  (X1+. 1 , 19 . 5 ) \  GRAPHIC  PRINT  'Log  S'\  RETURN 

880  SET  POSITION  (6,2.S>\  GRAPHIC  PRINT  USING  'EF-Ferai=  #♦#.###' ,EF 

881  SET  POSITION  <6,2.45>\  GRAPHIC  PRINT  USING  'Scrn  B  =  ###.##*', B 

882  SET  POSITION  (6,2.1 )\  GRAPHIC  PRINT  USING  'ED-Scrn  =  ###.### ',ES 

883  SET  POSITION  (6,1.75>\  GRAPHIC  PRINT  USING  'EV-GAP  =  ##*.#*#', EG 

884  SET  POSITION  (-.2,3)\  GRAPHIC  PRINT  3\  RETURN 
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